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ABSTRACT 

We have performed a statistical analysis of a sample of 28 nearby galaxy groups derived 
primarily from the Two-Dimensional XMM-Newton Group Survey (2dXGS), in order to as- 
certain what factors drive the observed differences in group properties. We specifically focus 
on entropy and the role of feedback, and divide the sample into cool core (CC) and non cool 
core (NCC) systems. This is the first time the latter have been studied in detail in the group 
regime. We find the coolest groups to have steeper entropy profiles than the warmest systems, 
and find NCC groups to have higher central entropy and to exhibit more scatter than their 
CC counterparts. We investigate the entropy distribution of the gas in each system, and com- 
pare this to the expected theoretical distribution under the condition that non-gravitational 
processes are ignored. In all cases, the observed maximum entropy far exceeds that expected 
theoretically, and simple models for modifications of the theoretical entropy distribution per- 
form poorly. A model which applies initial pre-heating through an entropy shift to match the 
high entropy behaviour of the observed profile, followed by radiative cooling, generally fails 
to match the low entropy behaviour, and only performs well when the difference between the 
maximum entropy of the observed and theoretical distributions is small. Successful feedback 
models need to work differentially to increase the entropy range in the gas, and we suggest 
two basic possibilities. 

We analyse the effects of feedback on the entropy distribution, finding systems with a 
high measure of 'feedback impact' to typically reach higher entropy than their low feedback 
counterparts. The abundance profiles of high and low feedback systems are comparable over 
the majority of the radial range, but the high feedback systems show significantly lower central 
metallicities compared to the low feedback systems. If low entropy, metal-rich gas has been 
boosted to large entropy in the high feedback systems, it must now reside outside 0.5 r5oo, 
to remain undetected in our study. Considering gas as a function of scaled entropy, we find 
similar levels of enrichment in both high and low feedback systems, and argue that the lack 
of extra metals in the highest feedback systems points to an AGN origin for the bulk of the 
feedback, probably acting within precursor structures. 
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1 INTRODUCTION 

The hierarchical view of structure formation indicates that present 
day clusters were assembled from galaxy groups. Indeed, th e ma- 
jority of galaxies in the Universe are situated in groups ( iTullvl 
fl987h . so understanding the processes at work in these systems is 
crucial to our understanding of galaxy properties as a whole. Inter- 
estingly, groups also seem to show a very diverse range of prop- 

* E-mail: ria@star.sr.bham.ac.uk 



erties (e.g. lOsmond & Ponmarill2004h . and understanding the root 
causes of these differences is key to understanding the processes at 
work in shaping present day galaxy clusters. It has also been known 
for some time that groups do not obey self-similarity, and for ex- 
amp le, show a steepening in the Lx — Tx relation compared to clus- 
ters jHelsdon & Ponma 320001) . This is explained via the increased 
importance of baryonic physics; cooling, feedback from AGN and 
supemovae, merger shocks and galaxy winds can have significant 
effects on the properties of the gas within the shallower potential 
wells of low mass clusters. 
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Knowledge of the dark matter-dominated potential, coupled 
with the entropy distribution of the gas within a system in hydro- 
static equilibrium, completely defines the therrnodyn amic proper- 
ties of the intracluster medium (ICM) dVoit et alj20"03h . The advent 
of high quality X-ray data from Chandra and XMM-Newton have 
allowed detailed studies of the entropy properties of clusters (e.g. 
IPratt. Amaud & PointecouteaJlioO^ : iMorandi & EttorillioOTh and 
simulations have been extensively used to predict properties such 
as entropy pro files and the chemical enrichment of the ICM (see 
iBorgani et al.l 2008a,b for recent reviews). However, the impor- 
tant group regime has been less well explored, although significant 
progress is now being made in modelling feedback effects on the 
group scale JPave. O ppenhei mer & Sivanandamll200 8). There is a 
parallel shortage of detailed obsenational evidence which bears on 
the processes shaping group evolution, and the aim of the present 
study is in part to address this. 

Observations of galaxy clusters show a striking bimodal- 
ity in their observed ICM properties, which leads to their clas- 
sification into cool core (CC) and non cool core (NCC) sys- 
tems. These classes are characterised by the observation (or 
not) of central positive temperature gradients. The open ques- 
tion in this field is how this apparent bimodality is produced. 
The short central cooling times observe d in both types of system 
jSanderson. Ponman & O'Sullivanl2006h . causes difficulties for at- 
tempts to explain the two classes within the same framework of 
cluster e volution. However, rec ent attempts have been made to do 
just this jMcCarthv et alj2008h . by invoking different levels of pre- 
heating before cluster collapse. 

Recent results show strong similarities between CC clusters 
and CC groups . For example, the ob s erved abundance gradients 
in CC clusters jPe Grandi et alj|2004l ; iBaldi et alj|2007h are also 
seen in the Chandra sample of 15 groups of Rasmussen & PormianI 
( l2007t) . 14 of which were observed to have CCs. However, one cru- 
cial area for understanding the role of feedback that has not yet 
been probed is the nature of NCC systems in the group regime. 

We are now entering an era where large studies of galaxy 
groups are p o ssible with high quality X-ray data jpinoguenov et al.l 



groups are p o ssible witn nign quality A-ray data arinoguenov et al 
120061 l20oi iRasmussen & Ponmaij |2007| ; ISun et alj |2008[) . Al 



though the sample analysed here was not selected in a statistical 
manner, we have an opportunity to gain considerable insight into 
the diversity of group properties, as this is the largest sample of 
groups to date with high quality XMM-Newton data, which has 
been analysed in a homogeneous way. One key advantage of our 
approach is that in the spectral analysis of this data, spherical sym- 
metry was not an a priori requirement (see lFinoguenov et al.l2006l 
l2007h . Additionally, due to the size of the sample, we have been 
able, for the first time, to separate the sample into CC and NCC 
groups, and also to look at the behaviour of groups as a function of 
temperature. This provides a unique opportunity to examine which 
processes drive the observed differences in group properties, by 
considering the behaviour which diverges from the mean relations. 

The layout of the paper is as follows. In Section[2]we describe 
the sample, analysis and basic properties of the groups, in Section|3] 
we divide the sample to look at the mean properties of sub-samples 
of the data, and Section|4]explores the radial properties of the sam- 
ple. We explore the entropy distributions and discuss the effects of 
feedback on the groups in Section [5] and we summarise our main 
results in Section|6] 



2 GROUP SAMPLE 

The Two-Dimensional XMM-Newton Group Survey (2dXGS) is 
an archival study of nearby (z < 0.024) galaxy groups, which 
were selected from the group catalogue of Mulchae v et alj ( 1200 3l) 
and were chosen to have pu blicly available XMM-Newton data 
dPinoguenov et al.l2006ll2007h . The complete sample selecti on and 
data analysis is described by Finoguenov et al.l ( l2006ll2007h . here- 
after referred to as F-l-06 and F-l-07 respectively. The former con- 
centrates on the 'low' redshift sample (z < 0.012), and the latter 
concentrates on the 'high' redshift sample (0.012 < z < 0.024). 
These works studied the radial properties of individual groups and 
their deviations from the mean profiles, applying a novel spectral 
approach which we briefly describe below. Our aim is to bring these 
two samples together , supplemented by groups from the sample of 
iMahdavi etal] ( l2005h which were analysed using the same proce- 
dures, in an effort to undertake a full statistical study of the prop- 
erties of a large sample of galaxy groups derived from high quality 
X-ray observations. 

The full 2dXGS sample contains 25 nearby galaxy groups ob- 
served with XMM-Newton, which have been analysed in a homo- 
geneous way. A summary of the data analysis procedure is given 
here; we refer the reader to F-l-06 and F-l-07, and references therein, 
for a full description of the XMM-Newton data, reduction and anal- 
ysis. There are two main stages in the analysis following the ini- 
tial XMM-Newton data reduction. Firstly, temperature and surface 
brightness maps were used to look at the overall structure of the 
group. The second part of the analysis extracted spectra from re- 
gions of contiguous surface brightness and temperature, and fitted 
single temperature hot plasma (APEC) models to the spectra, yield- 
ing the spectral properties of the group with no a priori assumption 
of the syste m being spherically symme tric. The abundances used 
were those o f lAnders & Grevessd ( fT98^ . and absorption was fixed 
at the Galactic value. 

The result of this approach is to yield a series of spectrally de- 
rived parameters in both two-dimensional regions, and also from a 
more traditional spectral analysis of a series of concentric annuli. 
We can look at the radial properties of each group by assuming the 
characteristic radius for the measurement to be the mean radius of 
the region. Here we will concentrate mainly on using the proper- 
ties derived from this novel analysis of 2d regions, except where 
otherwise stated. Due to the nature of the analysis, the customary 
spherical deprojection of the spectra cannot be performed. Instead, 
F-l-06 and F-l-07 determined the three dimensional gas properties in 
each of the analysed regions, by estimating the projection length of 
each region. We use here the derived gas properties from this ap- 
proach, and refer readers to F-l-06, F-l-07 and references therein for 
more details on this procedure. 

We define the mean temperature of each group to be that re- 
covered from fitting a spectrum extracted from the radial range 0. 1- 
0.3 rsoo (F-l-06, F-l-07), so we restrict the analysis presented here to 
the 21 2dXGS groups with measured temperatures in this range. We 
have supplemented the 21 groups in the 2dX GS sample wi t h seve n 
groups with the highest data quality from IMahdavi et al.l (l2005b, 
leading to a final sample of 28 galaxy groups. The Mahdavi et al.l 
(2003) groups were analysed using the same two-dimensional pro- 
cedure, but they have not also been analysed with the traditional 
annular approach which was also applied in the case of the 2dXGS 
groups. The groups in the Mahdavi sample were initially selected 
from the ROSAT All-Sky Survey/Center f or As trophysics Loose 
Systems (RASSCALS) of IMahdavi et^ ( l2000h . and cover red- 
shifts between 0.016 and 0.037. The mean temperatures for these 
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groups have been re-extracted, to cover the 0.1-0.3 rsoo region ap- 
phed to the 2dXGS groups. 

2.1 Group properties 

The basic properties of the group sample appear in Tables [T] and (2] 
the first of which refers to the 2dXGS sample of F+06 and F+07, 
and the seco nd of which shows the properties of the supplementary 
groups from iMahdavi et all j2005l) . Co lumn (1) shows the g roup 
name as expressed by F+06, F+07 and lMahdavi et alj ilOOS) . We 
now describe the origin of the values shown in Tables [T] and |2] We 
assume a Hubble constant of Ho = 70 kms^^Mpc^^ throughout. 

2.1.1 2dXGS groups 

Due to the overlap of 15 of the 2dXGS groups with the 
Group Evolution Multiw avelength Study (GEMS) sample of 
lOsmond & PonmanI j2004h . the latter is our primary source for 
the group properties (shown in Table [TJ supplementary to those 
provided by F +06 and F+07. Wher e distan ce measurements were 
available from lOsmond & PonmanI ( l2004h , these were used. For 
the remaining groups, we re-scaled the distance measurements of 
iMulchaev et alj j2003l) to our assumed Hubble constant (3C 499, 
NGC 507, NGC 2300 and NGC 4168). In the case of Hickson 51 
and the Pavo group, it was necessary to estimate the distance from 
the redshift information provided in F+07. 

X-ray lumino si ties are also drawn primarily from 
lOsmond & PonmanI j2004l) . which employed ROSAT PSPC 
data, and has the advantage of extending to larger radii than 
XMM-Newton data for many of these groups. We used the s u pplied 
/3-model surface brightness fits of lOsmond & PonmanI ^20041) 
to re-scale the X-ray luminosities from their original extraction 
ra dius to our rspp valu e s. Ind ividual extraction radii can be found 
in lOsmond & PonmanI ( |2004|) ; they are typically ~55 percent of 
our rsoo. In some cases, Osmond & Ponman could not fit a model 
to the surface brightness profile, and in these cases, a standard 
model of P = 0.5 and icorc = 6 kpc was assumed. This affects 
the following groups in our sample: Hickson 15, Hickson 92 and 
NGC 5171. F or groups not in the GE MS sample. X-ray luminosi- 
ties are from iMulchaev et alj ( 1200.^ . but have been corrected in 
the following way for differences in the assumed Hubble constant, 
and to extrapolate to rsoo in line with the X-ray luminosities 
of the GEMS groups . Where /3-mod e l surf ace brightness fits 
were available from IMulchaev et"ai] ( |2003|) (NGC 507), we 
re-scaled the X-ray luminosities both for the assumption here of 
a lower Hubble constant, and to scale from the original extraction 
radius of 200 hj^QQ kpc to our values of rsoo . In two cases (3C 449 , 
NGC 2300), the fitted core radii reported bv lMulchaev et al.l j2003h 
are lower than the resolution limit of the ROSAT PSPC instrument. 
In this case we set the core radius at this limit when scaling the 
X-ray luminosity to our value of rsoo. For NGC 4168, w here only 
an upper limit on Lx was determined by IMulchaev et al.l j2003l) 
and hence no surface brightness modelling was undertaken, we 
report the upper limit corrected for the different Hubble constant 
only. In the case of Hickson 51, which is not in the GEMS sample, 
we used the X-ray lu minosity for the group NRGbl28 from the 
RASSCALS sample of lMahdavi et"aLl ( l200(i) . which is given as an 
equivalent designation for Hickson 51 in the NASA/IPAC Extra- 
galactic Database (NEDfl We extrapolated this X-ray luminosity 
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from the original extraction radi us of 0.5 hTnn Mpc t o our rsoo 
assuming the standard /3-model of lOsmond & PonmanI ^20041) . and 
we also scaled the X-ray luminosities to a Hubble constant of Hq 
= 70 kms^^Mpc^^. This is the same procedure w e have used to 
determ ine X-ray luminosities for the majority of the lMahdavi et al.l 
( l2005l) groups (see Section |2. 1.21 ). so there is a consistency in our 
approach. 

It was necessar y to estimate the X-ray luminosity of Pavo 
from the work of Ma chacek et "d]( l2005h . who fitted a (3-mode\ sur- 
face brightness profile to the extended emission outside the central 
galaxy (NGC 6876) with r^ore = 196" and (3 = 0.3, yielding an X- 
ray luminosity (in the energy band 0.5-2 keV) of log Lx = 41.8 
within 463", assuming a Hubble constant of 75 kms~^Mpc~^. We 
have extrapolated this model to our rsoo and have corrected for the 
differences in the assumed Hubble constant, to quote the X-ray lu- 
minosity for Pavo shown in Table[T] The remaining caveat with this 
value is the slightly different energy band compared to the remain- 
der of the sample, which were all derived from ROSAT data. 

The velocity dispersion of each group comes from F+06 and 
F+07, except in the case of NGC 4636, for which F+06 could not 
determine the group membership satisfactorily (we refer the reader 
to this work for more information), and in this c ase, we quote the 
velocity dispersion of iOsmond & PonmanI ( |2004|) . The mean tem- 
peratures of the groups and the values of rsoo were all derived by 
F+06 and F+07 using the method described in Section|2| 

We also show selected optical properties of the central galaxy 
and the group. Values of D25, the diameter of the isophote where 
the B-band surface brightness is 25 mag arcsec"^, are deter- 
mined for the central group galaxy from the RC3 catalogue of 
4e Vaucouleurs et al. (1991). We have applied the procedure of 
lOsmond & PonmanI (12004) in determining the B-band luminosity 
of the brightest group galaxy (BGG) and of the group, by extract- 
ing galaxies from NED within a projected ra dius of rsoo, centred 
on the group co-ordinates. For groups in the lOsmond & PonmanI 
(2004) sample, we used the group co-ordinates provided in this 
work, and for the remaining groups we use the NED co-ordinates 
corresponding to the group name. The group luminosities are 90 
percent complete and we applied an absolute magnitude cut of Mb 
= -16.32 to the group galaxies. The BGG is chosen to be the bright- 
es t galaxy within 0.25 rsoo of the group centre. We refer the reader 
to lOsmond & PonmanI ( |2004|) for more information on the applied 
method. In the ca se of the Pavo group, w e have assumed the BGG 
to be NGC 6876 dMachacek et alj|2005l) : this is confirmed by our 
procedure for extracting B-band luminosities. Table [T] also gives 
the maximum radius in units of rsoo to which spectral information 
is available, for each of the groups, to indicate the completeness 
of the spectral coverage. The final columns in Table [T] denote the 
subsamples to which each group belongs, in terms of their mean 
temperatures and core properties. These classifications are defined 
in Section|3] 

2.1.2 \Mahdavret aZI (l2005l) groups 

Three of the Mahdavi et aP ( I2OO5I) groups appear in the sample of 
lOsmond & Ponman C2(}04). Thes e are Hickson 97, SRGbl 19 (cen- 
tral galaxy NGC 741 : iMahdavF et al. 200 ^ and NGC 5129. For 
these groups, distances are from lOsmond & PonmanI ( l2004h . but 
in the remainder of c ases the distances we re estimated from the 
redshifts presented bv iMahdavi et ai] ( I2OO5I) . Similarly, the X-ray 
luminosities for these 3 grou ps are re-scaled from the or iginal ex- 
traction regions presented in lOsmond & Ponmad ( I2OO4I) to match 
our rsoo values, as described in Section l2.1.1l for the groups in the 
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Table 1. The basic properties of the 2dXGS grou ps in the sample. The groups have been classified by their mean temperature and their core properties (see 
Section |3)- Distances and values of Lx are from lOsmond & PonmanI (2004) unless otherwise stated. For data from Osmond & Ponman (2004), the X-ray 
luminosities have been extrapolated to our rsoo values (see text). Lg bqq and Lg g^p are calculated as described in the text. Mean temperatures and 
values of rgoo are from F-l-06 and F+07. Velocity dispersions are also from F+06, F-l-07, except in the case of NGC 4636, where the value used is that from 
[Osmond & Ponman (2004), and values of D25 (the diameter of the isophote where the surface brightness is 25 mag/arcsec^ in the B-band) are from RC3 
tde Vaucouleurs et al.i,1991]) , and refer to the brightest group galaxy (BGG). Vmax denotes the maximum radius to which spectral information is available, in 
units of rsoo . 



Group Dist. log Lx a T rsoo f 25 log Lb log Lb r-max Warm/Cool CC/NCC 

(Mpc) (ergs^i) (kms^^) (keV) (kpc) (arcmins) {1^b,q) (Ls,©) {rsoo) 

BGG Group 





D / 


AO Q'2_l_0.26a 

42.93±o 27 


^^^4-1 1 
JJJitiiZ 


1 oQ_i_n m 

1 .ZoitU.UZ 


4J J 


1 1 n 


1 A 7 A 


1 A 77 


A A'\ 


Warm 






yj 


4Z. IVitU.Uj 


/I A/1 _L 1 TO 


U.oZitU.U4 


Zoo 


U.a/ 


1 n A/1 
1U.U4 


1 A 


A <7 
U.jZ 


L-OOl 






04 


/IT n<_i_n m 


Z(5Zlt4J 


U. / JltU. Ly 


jZ4 


Z.yj 


1 A 0^ 
iU.VJ 


1 1 77 
i i .Z / 


A 7S 


1..001 




HCG 51 


114* 


'~T^.\J^ 1_ \J . 1 _J 


546±151 


1.16it0.13 


424 


1.10 


10.51 


11.10 


0.82 


Warm 


NCC 


HCG 62 


74 


43.18±0.04 


418±51 


1.06±0.02 


403 


_d 


10.54 


11.32 


0.46 


Warm 


CC 


HCG 68 


41 


41.71±0.04 


191±68 


0.69±0.09 


308 


2.19 


10.63 


11.32 


0.19 


Cool 


NCC 


HCG 92 


88 


42.16±0.04 


467±176 


0.79±0.24 


334 


1.91 


10.51 


11.01 


0.42 


Cool 


NCC 


IC 1459 


26 


41.42±0.04 


256±|| 


0.59±0.03 


280 


5.25 


10.63 


10.90 


0.21 


Cool 


NCC 


NGC 507 


70" 


43.37±0.03" 


580±94 


1.34±0.01 


467 


3.09 


11.00 


11.60 


0.53 


Warm 


CC 


NGC 533 


76 


42.71±0.03 


439±60 


1.26±0.01 


448 


3.80 


10.99 


11.39 


0.56 


Warm 


CC 


NGC 2300 


29" 


41.93±0.04" 


278±|f 


0.75±0.01 


339 


2.81 


10.34 


10.67 


0.32 


Cool 


CC 


NGC 2563 


73 


42.58±0.03 


384±49 


1.3 1 ±0.05 


460 


2.09 


10.64 


11.35 


0.36 


Warm 


CC 


NGC 4073 


96 


43.46±0.02 


565±72 


1.87±0.05 


575 


3.16 


10.96 


11.65 


0.49 


Warm 


CC 


NGC 4168 


38" 


<40.87" 


259±|| 


0.77±0.31 


258 


2.75 


10.39 


10.81 


0.19 


Cool 


NCC 


NGC 4261 


41 


42.33±0.03 


429±i4 


1.11 ±0.02 


451 


4.07 


10.85 


11.46 


0.19 


Warm 


CC 


NGC 4325 


117 


43.16±0.01 


376±70 


1.01±0.01 


389 


0.95 


10.65 


10.97 


0.54 


Cool 


CC 


NGC 4636 


10 


41.72±0.02 


284±73 


0.77±0.01 


331 


6.03 


10.02 


10.39 


0.18 


Cool 


CC 


NGC 5044 


33 


43.10±0.01 


357±|« 


1.21±0.01 


430 


2.95 


10.50 


11.15 


0.40 


Warm 


CC 


NGC 5171 


107 


42.48±0.06 


494±99 


1.21±0.05 


436 


1.10 


10.77 


11.44 


0.49 


Warm 


NCC 


NGC 5846 


30 


42.01 ±0.02 


368±5i 


0.69±0.01 


309 


4.07 


10.73 


11.10 


0.13 


Cool 


CC 


Pavo 


57'' 


42.69*= 


440±96 


0.77±0.12 


330 


2.82-f 


11.02 


11.47 


0.49 


Cool 


NCC 



Notes: 



" Values from lMulchaev et aP i2003h . rescaled using Ho = 70 kms~^ Mpc~^ and extrapolated to our rsoo values (see text for details). 

Distance estimated from redshift given in lFinoguenov et al] i2007l) . 
" Value from [ M ahdavi e t alj j2000l) rescaled using Hq = 70 kms~^ Mpc~^ and extrapolated to our rsoo values using the standard /3-model of 
lOsmond & Ponmarj j2o64r (see text for details). 
''■ Not available in RC3. 

" Lx extrapolated to rsoo and corrected for Hq = 70 kms~^ Mpc^^ from the measured X-ray luminosity and suiface brightness profile of lMachacek et al] 
i2005l) (Note: energy band is 0.5-2 keV for this value). 

■'' D2S from RC3 assuming central galaxy is NGC 6876 jMachacek et al.l2005l) . 



2dXGS sample. For the remaini ng four groups, we use the X-ray lu- 
minosities from lMahdavi et al] (12000,) . which were extracted from 
the ROSAT All-Sky Survey, with in a radius of 0.5 hr o pMpc . As- 
suming the standard /3-model of lOsmond & Ponma we 
scale these X-ray luminosities both for Ho = 70 kms~^Mpc^^ 
and to our radius of rsoo, in the manner described in Section 
|2T7] This ensures a dec ree of consistency between the 2dXGS 
and lMahdavi et all ([2OOS) samples. 

Velocity disper sions for all seven systems come from 
iMahdavi et alj J2005j), and mean temper atures and rsoo values 
were re-derived from lMahdavi et al ]l l2005h in the radial range 0. 1- 
0.3 rsoo- The B-band optical luminosities of the BGG (Lb.bgg) 
and the group (Lb, group) were determined using the same proce- 
dure as in Section l2.1.1l For the three GEMS groups , we used the 
co-ordinates provided bv lOsmond & Ponm anI (l200 4i) for th e NED 
search; we use the group co-ordinates from Mahda vi et alj j2005l) 
for the remainder of the sample. iMahdavi et al. ( 2005) indicate that 
Abell 194 has no dominant galaxy, but our method finds NGC 541 
to be the brightest galaxy within 0.25 rsoo, and we a dopt this as the 
BGG . The central galaxy of RGH 80 (NGC 5098: IMahdavi et al.l 
I2OO5I) is listed in NED as a galaxy pair, but we find PGC 046515 to 
be the brightest galaxy within 0.25 rsoo of the group centre, and we 



adopt this galaxy as the BGG in this system. Values of D25 were 
again extracted from the RC3 catalogue of Ide Vaucouleurs et al.l 
( I199 ll) for these systems. Table |2] also gives the maximum radius 
in units of rsoo to which spectral information is available, for each 
of the groups, to indicate the completeness of the spectral coverage. 



3 DIVIDING THE SAMPLE 

The 2dXGS sample of groups has been presented previously (F-l-06, 
F-l-07) in two subsamples, distinguished by group distance, but it 
is also instructive to partition the sample on the basis of physical 
properties. We have divided the groups into two main subsamples, 
firstly on the basis of the mean temperature of the group, which acts 
as a proxy for the group mass, and secondly on the basis of whether 
or not the group has a CC. The classification of galaxy clusters into 
CC and NCC systems is well-known (e.g. Peres et al. 1998). These 
systems have been shown to exh ibit different observational proper- 
ties (e.g. ISanderson et alj|2006h , and henc e it is instructive to ex- 
tend th is classification to the group regime. iRasmussen & PonmanI 
1 I2OO7I) studied the temperature and abundance profiles of 15 groups 
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Tabl e 2. The basic p roperties of the groups in the sample originally from lMahdavi et al] j2005l) . Mean temperatures and values of rsoo have been re-derived 
from iMahdavi et all (2005) in the radial range 0. 1-0.3 rsoo- Di stances have been estim ated from th e redshifts given in Mahdavi et al. ( 2005), and X-ray 
luminosities have been extrapolated to our rsoo values from Osmond & Ponman 1 2004) and iMahdavi et aU (2000). Values of Lg bqq and L^^grp ^u'e 



calculated as described in the text, 
is available, in units of rsoo . 



D25 values are from RC3 i de Vaucouleurs et al. 1991). Vmax denotes the maximum radius to which spectral information 



Group 


Dist. 
(Mpc) 


log Lx 
(ergs-l) 


A 194 


76" 


42.88i0.05'' 


HCG 97 


92 


42.77±0.10 


NGC5129 


108 


43.14±0.09 


NRGb 184 


96" 


42.57±0.12'' 


RGH80 


158" 


43.18±0.06'' 


SRGb 119 


79 


42.65±0.10 


SS2bl53 


68" 


42.69±0.07'' 



(kms-l) 



T 


''"500 


£'25 


log Lb 


log Lb 


f'max 


CoolAVarm 


CC/NCC 


(keV) 


(kpc) 


(arcmins) 


(Lb.o) 
BGG 


(Lb.q) 
Group 


(rsoo) 






1.01±0.15 


393 


1.78<= 


10.74': 


11.56 


0.86 


Cool 


NCC 


1.20±0.05 


439 


1.66 


10.39 


11.05 


0.69 


Warm 


CC 


0.95±0.03 


379 


1.70 


11.05 


11.37 


0.53 


Cool 


CC 


1.37±0.09 


477 


0.87 


10.41 


11.04 


0.68 


Warm 


CC 


1.16±0.02 


429 


0.78'* 


10.59'' 


11.37 


1.4 


Warm 


CC 


1.34±0.07 


470 


2.95 


11.11 


11.29 


0.58 


Warm 


CC 


0.83±0.01 


348 


2.09 


10.61 


10.89 


0.68 


Cool 


CC 



550±«« 
383±50 
283±2f 
390±|0 
602±|f 
416±i| 
161±i? 



Notes: 

Distances estimated from redshifts given in IMahdavi et aP i2005l) . 
' Values from Mahdavi et al. 12000) extracted within 0.5 Mpc, rescaled here for Ho = 70 kms~^ Mpc~^ and extrapolated to our rsoo using the standard 
/3-model of Osmond & Ponman 1 2004) (see text for details). 

NGC 541 assumed as BGG (see text for details) as no dominant galaxy Hsted bv lMahdavi et aljj2005l) . 
'' PGC 0465 15 assumed as BGG (see text for details). 



with Chandra data, 14 of which were noted to have a CC. The prop- 
erties of groups that do not have CCs, and the connection between 
the overall state of the ICM and the core properties of groups, has 
not been investigated to date. A sample of the size presented here, 
although not statistically selected, offers an opportunity to probe 
the wider properties of groups compared to the presence (or not) of 
aCC. 

To enable us to compare the radial profiles of groups directly, 
we scale the radii by rsoo, the radius within which the mean density 
of the group is 500 times the critical density. This was calculated in 
an iterative fashion using (F-l-06, F-l-07), 

rsoo = 0.39lf "■''^/ifoV (1) 

The partition of the groups into subsamples based on their mean 
temperature and core properties is detailed below. 

3.1 Temperature 

We define the mean temperature of each group as that determined 
within the radial range 0.1-0.3 rsoo (F-l-06, F-l-07). This removes 
the effect of any CC on the temperature determination, providing a 
robust measure of the mean temperature of the system. The median 
temperature across all groups was found to be 1.035 keV, and this 
was used to divide the groups into 2 subsamples, the Cool groups 
with f < 1.035 keV, and the Warm groups with f > 1.035 keV. 
The temperature classification for each group appears in Tables [Tj 
andU 

3.2 Core properties 

To distinguish between CC and NCC systems, one tra ditional ap- 
proac h is to cons ider the cooling time of the system (e.g. lPeres et al.l 
Il998h . However, I Sanderson et alj J2OO6I ) showed that in the cluster 
regime, systems with no central cooling can also exhibit the short 
central cooling times (< 5 Gyr) typically associated with CC sys- 
tems. They defined the presence of a CC by considering the ratio 
of temperatures between an inner radius and an intermediate radial 
range, and we adopt the same approach here. 

We define a group to have a CC if the mean temperature in 



the region 0.1-0.3 rsoo is greater than the mean temperature in the 
range 0.00-0.05 rsoo, indicating a positive temperature gradient in 
the inner regions. If this was not the case, the group is classified 
as a NCC system. This classification is based wherever possible 
on results from the analysis of both the annular and the 2d regions 
reported in the earlier 2dXGS papers, to improve the radial cover- 
age and allow a more thorough appraisal of the temperature profile. 
For the seven groups from the Mahdavi et al. (2005) sample, only 
the 2d information is available. The temperature profiles of the CC 
systems are shown in Figure [Tj and those for the NCC systems are 
shown in Figure|2| In both cases, a weighted local regression fit to 
the data is also shown, to indicate the overall behaviour of the tem- 
perature profile whilst suppressing scatter. This fit was performed 
using the 'LOESS' fu nction in version 2.5.1 of the R sta tistical en- 
vironment packag^f] IR Development Core Teadl2008 f), hereafter 
referred to as R. The LOESS algorithm performs a weighted least 
squares fit in the local neighbourhood of each data point, where the 
size of the neighbourhood is defined to include a specified propor- 
tion of the data sample, and the distance to each neighbour is used 
to weight the least squares fit. For more information on the LOESS 
algorithm, we refer the reader to Cleveland, Grosse & Shyu ( 19921). 

For two systems, classification on the basis of the above cri- 
terion was problematical. Hickson 92 would be designated a CC 
system, since its profile shows a small drop in central temperature, 
albeit of rather low significance. However, this group (commonly 
known as Stephan's Quintet) is a system of galaxies undergoing 
multiple collisions, w hich are currently disturbing and heating its 
intergalactic medium dTrinchieri erZll2003l I2OO5I) . It is far from 
being a typical CC system, in which radiatively cooling, dense gas 
is centred on a dominant early-type galaxy. We therefore reclas- 
sify it as an NCC system. In the case of Abell 194, the quality 
and resolution of the data did not permit us to extract a central 
spectrum within 0.05 rsoo- However, the temperature profile rises 
consistently inward to the innermost point, at ~0.06rsoo- We are 
therefore confident in classifying this group as a NCC, since all our 
CC systems show CC behaviour by this radius, whereas Abell 194 
does not. The final sample therefore consists of 10 NCC groups, 

^ http://www.r-project.org 
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Figure 1. The temperature profiles of the groups classified as being CCs. In the case of the 2dXGS groups, the data from both the annuli and two-dimensional 
regions has been used. The solid lines show weighted local regression fits to the data, see the text for details. The plots are arranged in order of ascending mean 
temperature, from top left to bottom right. Red tick marks show the radius of maximum temperature used in the calculation of the temperature drop, and the 
dashed red line shows a straight line fit (fitted in log-log space) to the temperature profiles in the core region. Continued overleaf. 



and 18 CC groups. Ten of the l atter a re included in the Chandra 
study of lRasmussen & PonmanI i2007h . 

One interesting property of CC groups is the magnitude of the 
drop in temperature within the core. To estimate this, we found the 
radius of peak temperature, from the maximum in the local regres- 
sion fits to the data. To be able to evaluate the temperature drop to 
the same small radius (0.01 rsoo) in all groups, we fitted a straight 
line in log-log space to the temperature within the radius of max- 
imum temperature (red dashed lines in Figure [T), forcing the con- 
dition of passing through this maximum point. This fit was then 
used to evaluate the temperature at 0.01 rsoo, allowing groups with 
different radial sampling to be directly compared. 

For our CC subsample, the ratio of the temperature at 
0.01 rsoo to the mean temperature of the system is found to be 



Tg/T = 0.69±0.16. This ag rees with the result (0.58±0.14) of 
iRasmussen & PormianI 1 20071) within errors. Ho wever, the 'central 



temperatures' of ( IRasmussen & PonmanI l2007h were measured in 
the innermost radial bin of their analysis, which in most cases lay 
within 0.01 rsoo - This may account for their slightly lower value of 

We quantified the observed temperature drop as the differ- 
ence between the maximum temperature and that interpolated at 
0.01 rsoo - This quantity is temperature dependent, as hotter systems 
have the capacity to show larger temperature drops. To remove this 
dependence, we divided the temperature drop by the peak tempera- 
ture of the system. Panel (a) of Figure [3] shows this fractional tem- 
perature drop versus mean temperature. We can test for possible 
correlation here using Kendall's rank order correlation coefficient 
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r, which is found to be 0.15 with a p- value (probability of chance 
occurrence) of 0.38, indicating no significant correlation within our 
sample. However, this is partly due to the low value of the central 
temperature decline in our hottest system, NGC 4073, and exclud- 
ing this system yields r = 0.27 and a p- value of 0. 14, indicating a 
positive correlation at the >lcr level. The values of the fractional 
dip in Figure[3ja) are clearly for the most part low er than the value 



0.6 w hich is typical in CC clusters (see F ig. 7 of [ Sanderson et al 
l2006t) . in agreement with the findings of iRasmussen & Poiunan 
( 120071) . The CC group with the lowest fractional temperature drop 
(Tdrop I Tp eak '•^ 0.1) is SS2bl 53, the central galaxy of which is 
NGC 3411 jMahdavi et"^l2005l) . This group has been interpreted 
as evide nce of a CC s ystem that has been re-heated by recent AGN 

activity jO'SuUivan e t al. 2007J^ 

The Chandra study of IRasmussen & PormiarJ ('2007) found 
that the radius at which the temperature peaks is correlated with 
group temperature (and hence group size). We test for a correla- 
tion between the physical sizes of the CCs and the value of rsoo 
(see Figure [3] panel (b)) in our sample. The value of Kendall's 
rank order correlation coefficient r is 0.29, with a p-value of 
0.10, so we have some evidence for the existence of a correla- 
tion: larger systems tend to have larger CCs, which is perhaps not 
a sur prising result, and agrees qualitatively with the analysis of 
IRasm ussen & Ponman (2007). Fitting a linear model using an or- 
thogonal regression dlsobe et alil990 ) to allow for intrinsic scatter, 
we find a relationship between the size of the CC and the value of 
''500 of the following form. 



r(r™,) = (0.14 ± 0.07)r5oa - (1.0 ± 15.4) kpc. 



(2) 



where r{Tmax) is the radius of maximum temperature in kpc, and 
rsoo is also m easured in kpc. This rela t ion yi elds a flatter slope than 
the analysis of lRasmussen & PonmanI ( l2007h . who found a slope of 
0.20±0.02 rsoo, although the large error bars mean that the slopes 
of the two relations agree within errors. There are some differences 
in the methodology used to estim ate rsoo in the two studies, with 
theH asmussen &PonmaiJ ( f2007l) values for the groups in common 
tending to be larger, which could contribute to the slightly different 
results obtained. 



Our result has essentially zero intercept, so that 
r{Tmax) ~ 0.14 rsoo. This is very similar to the ratio of the tem- 

S ieratu re peak radius to rsoo found for clusters bv I Sanderson et af] 
200g) . Hence it seems that although the depth of the CC is less in 
groups than clusters, their size (scaled to rsoo) is similar. In Figure 
[3] panel (c), we scale the temperature peak radius to rsoo and find 
no significant temperature trend in the scaled size of the CCs. 

Comparing the global properties of the CC and NCC groups 
in our sample, the mean temperature of the NCC systems is found 
to be 0.89±0.08 keV, compared to l.ll±0.07keV for the CC 
groups. Therefore, although the two populations overlap, the cen- 
troid of the distribution is significantly lower for the NCC subsam- 
ple. Similarly, the mean X-ray luminosities of the subsamples are 
log Lx = 42.70±0.12 erg s^^ for the CC groups, and log Lx = 
42.20±0.21 erg s~^ for the NCC systems. Of course, we are not 
dealing with a statistical sample of groups, so it is important not 
to over-interpret these results. However, one might expect selec- 
tion effects to work in the opposite direction — cool, faint groups 
will be easier to detect if they have CCs — so our results provide 
tentative evidence that cooler groups are more likely to lack CCs. 
Moreover, it is interesting to note that three of the four coolest sys- 
tems, with temperatures less than 0.7 keV, do not have CCs, which 
hints at the possibility of a lower temperature threshold for CC sys- 
tems. Further study of low temperature groups (which are faint, and 
therefore hard to detect) is required to confirm whether there is in- 
deed such a threshold. 

To summarise, for our sample the fractional central drop 
in temperature is smaller than that seen in clusters. We see no 
clear trend in CC strength with mean temperature within our 
group sample, but the coolest systems (T < 0.7 keV) tend not 
to show CCs at all. We find a relation rather similar to that of 
IRasmussen & PonmanI ( |2007[) between the physical size of the CC 
region (defined as the radius of maximum temperature) and rsoo. 
Smaller groups (lower rsoo) exhibit smaller CCs, and the core size 
across our temperature range is typically 14% of rsoo, in good 
agreement with what is seen in clusters. 
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Figure 2. The temperature profiles of the groups classified as being NCCs. In the case of the 2dXGS groups, the data from both the annuli and two-dimensional 
regions have been used. The plots are arranged in order of ascending mean temperature, from top left to bottom right. The solid lines show weighted local 
regression fits to the data, see the text for details. We note here that no smooth fit could be found for the temperature profile of NGC 4168, hence a fit is not 
shown in this case. 



4 RADIAL PROFILES OF GAS PROPERTIES 

Here we look at the radial variation of deprojected gas properties: 
entropy, pressure and gas density. These quantities have been de- 
rived from a two dimensional analysis, rather than an annular ap- 
proach, so we define the characteristic radius of each measurement 
as the mean of the two bounding radii of the spectral extraction 
region. 

4.1 Entropy profiles 

The entropy of a relaxed system, coupled with knowledge of the 
potential well , completely defi nes the properties of the intracluster 
medium (e.g. IVoit et allllOOSl) . Gas entropy can provide useful in- 
sights into the non-gravitational processes shaping the ICM, since it 



remains unchanged when gas is simply moved around within a sys- 
tem. W e define entropy as (e.g. lPormian. Sanderson & Finoguenovl 
l2003h . 

S = ^n^'/^ (3) 

where T is the mean temperature of the system, and Ue is the gas 
density, and S has units of keV cm^. This is related to the true 
thermodynamic definition of entropy via a logarithm and additive 
constant, and has the benefit of acting as a proxy for entropy which 
follows directly from two X-ray derived properties. 

In the self-similar case, entropy simply scales with the virial 
temperature of the system. However, a modified entropy scaling of 
5* / T'^/^ was found to perform well across a wide range of virial 
masses, from groups to clusters jPonman et al.ll2003l) . We aim here 
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Figure 3. (a): The fractional temperature drop (Tdrop^'^peak) in CC 
systems versus the mean temperature of the system, (h): The radius of the 
maximum temperatur e versus rspp for the CC groups (open circles). The 
dashed line is the fit of lRasmussen & PonmanI |2P07) to their sample of CC 
groups (excluding NGC 4125) analysed using Chandra data, with the dot- 
ted lines showing the upper and lower bounds of the fit. The solid line is 
an orthogonal regression fit to our data (see text), with its 1 a associated 
confidence region. The solid square marks the mean of the rsoo values for 
the NCC groups, and the associated error bar marks the minimum and max- 
imum of the rsoo values for the NCC groups in the sample. ( c): The radius 
of maximum temperature (as a fraction of rsoo ) versus the mean temper- 
ature of the system. The solid square marks the mean of the rsoo values 
for the NCC groups, and the associated error bar marks the minimum and 
maximum of the r'soo values for the NCC groups in the sample. 



to test the use of this modified entropy scaling across the group 
regime. Considering first the self-similar scaling (S oc T), Figure|4] 
shows power-law fits (in log-log space) to these scaled entropy pro- 
files for both the Cool and Warm sub-samples, with their associated 
68% confidence regionflf]. We note that few groups contribute be- 
yond ~0.5r5oo, and hence results shown here beyond this radius 
should be treated with caution. If the self-similar scaling works 
then it should remove any systematic temperature dependence in 
the profiles, so that the fits to the Cool and Warm samples should 
coincide. In practice, one can see that outside the core the scaled 
profile for Cool groups lies significantly above that for Warm ones, 
indicating that a scaling less strong than S oc T is required. 



4.1.1 Optimum entropy scaling 

We can investigate the optimum entropy scaling by looking at the 
entropy at a particular radius as a function of temperature. To each 
individual (unsealed) entropy profile, we fitted a power law model, 
to determine the entropy at 0. 1 rsoo • In Figure|5]we plot the entropy 



Throughout this section, the estimation of confidence regions on power- 
law fits uses a modification of the 'confidence. band' function written for the 
R statistical environment package by Derek Young and David Hunter, see 
rhttp://www.stat.psu.edu/~dhunter/R/confidence.band.r| 



Figure 4. Power-law fits to the entropy profiles for all groups in the sam- 
ple, scaled self-similarly, as 5 / f . The white dotted Hne and associated 
68% confidence region show a power-law fit to the Cool groups, and the 
black dashed line and associated 68% confidence region show a power-law 
fit to the Warm groups. The grey data points show individual entropy mea- 
surements for the Cool groups (open cir cles) and the Warm gro ups (open 
squares). The theoretical slope of r^ l iTozzi & Normarj|200lh showing 
the expectation from gravitational processes only is also shown. 




1.5 2 3 4 
Mean Temperature (keV) 

Figure 5. Entropy at 0. 1 rsoo as a function of the mean temperature of 
the system (measured between 0.1 a nd O.Brsop) for t he 2d XGS sample 
(solid squares). The cluster sample of lSanderson et al] i2009h is shown as 
simple error bars. The soli d line shows a BCES orthogonal regression fit 
<Akritas & Bershadvl [1993) performed in log-log space to the combined 
group and cluster sample. 



© 2009 RAS, MNRAS OQO.fTtilil 



10 R. Johnson, T. J. Ponman and A. Finoguenov 



at this characteristic radius as a function of the mean temperature 
of the system. 

We fit a straight line in log-log space to yield the optimum 
entropy scaling, which will simply be the slope of this relation. 
Due to the measurement errors on the mean temperature, and the 
errors on the entropy measurements, which are propagated from 
the errors on the power-law fits to each individual group, we fit the 
relation using the BC ES generalisation of t h e ort hogonal regression 
method presented by lAkritas & Bershad\] ( Il996) . This allows for 
both measurement errors (which may be correlated) and also for 
unknown intrinsic scatter in the relation. The narrow temperature 
range covered by the groups, coupled with the large dispersion in 
group properties, results in a large error in the slope of the S-T 
re lation. We therefore co mbine our sample with the cluster data 
of ISanderson et alj ( I2OO9I) and fit the relation using the combined 
group and cluster sam ple, to increase the temperat ure b aseline. 

As described bv lAkritas & Bershadvl ( 1 19961) and llsobe et al.l 

( ll99Cl() . the analytic form presented for the estimates of variance 
on the slope and intercept are only suitable for large samples, and 
more appropriately a bootstrap error analysis should be applied 
to describe the error behaviour of small samples. We have calcu- 
lated both the analytical variance from Akritas & Bershadv ( 1996|), 
and the result from a bootstrap analysis. We extracted 100 boot- 
strap samples from our original data, and calculated the slope and 
intercept in each case, determining the mean and standard devia- 
tion from these results. Applying the orthogonal regression method 
yields a slope of 0.79±0.06 using the analytical variance; the slope 
resulting from our bootstrap analysis was found to be 0.78±0.06. 

The recov ered slopes are a little steeper than the 2/3 scaling of 
iPonman et al. I1I2OO3), although they are consistent with it at the 2 a 
level. Our result is in good agreement with recent work of^ Sun et"al] 
( I2OO8I) , who find a slope of 0.78±0. 12 at the somewhat larger ra- 
dius of 0.15 rsoo, when applying a BCES y on x regression to their 
group sample . Since the S oc t"^!^ modified entropy scaling of 
IPonman et al. I (§003) provides an acceptable representation of the 
en tropy scaling fro m both the present study, and the Chandra study 
of ISun et"ai] ( l2008h . we adopt this scaling for the remainder of this 
paper. However, we should bear in mind that the actual scaling now 
seems likely to be rather steeper than T^/'^. 



4.1.2 Cool and Warm samples 

Adopting the modified entropy scaling as discussed above, to 
remove the mean effects of group temperature on the entropy, 
we work with the scaled entropy S j 'f'^l'^ ^ which has units of 
keV^''^ cm'^. Power-law models are fitted to radial profiles of the 
scaled entropy for the Cool and Warm subsamples separately. The 
resulting power-law fits and their error envelopes are shown in 
Figure |6] Where the entropy is only driven by gravitational pro- 
cesses, its radial variati on is expected to be proportional to r^'^ 
jTozzi & Normaj |200lh . Figure [6] shows the observed profiles to 
be much flatter than this, as also found by Sun et al. (200l). This 
is consistent with observations of the inner regions of NCC clus- 
ters, although CC clust ers are found to lie closer to the r^'^ relation 
dSanderson et alj|2009l) . The slopes of the entropy profiles fitted to 
the Cool and Warm groups are 0.77±0.03 and 0.65±0.02 respec- 
tively, indicating a steeper slope in the case of the cooler systems. 
Figure[6]also shows the 68% confidence bounds on the fits, arising 
from errors in the slope and intercept. It can be seen that outside the 
core, the scaled entropy for the Cool systems lies somewhat above 
that for the Warm groups, confirming the discussion from the pre- 




Figure 6. Power-law fits to the scaled (5 / f ^/^) entropy profiles for all 
groups in the sample. The white dotted line and associated 68% confidence 
region shows a power-law fit to the Cool groups, and the black dashed line 
and associated 68% confidence region shows a power-law fit to the Warm 
groups. The grey data points show individual scaled entropy measurements 
for the Cool groups (open c ircles) and the Warm g roups (open squares). 
The theoretical slope of r^ l jTozzi & Normanl200lll expected from gravi- 
tational processes alone is also shown. 



vious section, that the entropy scaling is actually rather steeper than 

5" oc T^/^ 



4.1.3 Core properties 

We now consider the entropy profiles of the groups, divided into the 
CC and NCC samples d escribed in Section [5] For a sample of CC 
clusters, iDonahue et al.l ( I2OO6) find that the addition of a constant 
entropy pedestal of ~10 keV cm'^ to a power-law entropy profile 
provides a better fit to their data than a simple power-law model. We 
use the notation 5*' to refer to scaled entropy, i.e. S' is equivalent 
to 5 / "f^l'^ , The pedestal modification is of the following form, 

y(r) =5-^ + 5^00 f — )°, (4) 

where So is the entropy pedestal, S500 is the normalisation of the 
profile at rsoo, and a is the index of the fit. The entropy pedestal 
modification of Donahue et al.l ( I2OO6I) is apparent within ~10 kpc; 
the corresponding radial range in our data is ~ 0.02-0.04 rsoo. 
However, within this radius there does not appear to be a flattening 
of the entropy profiles. This is exhibited in Figure|7] which shows 
(unweighted) local regression fits to the entropy profiles of the CC 
and NCC groups. We employed the 'LOESS' algorithm in R to per- 
form the fit, and fitted the data in log-log space. The profiles of the 
CC groups continually decrease inwards towards a scaled entropy 
of ~ 8 keV^^^ cm^ at radii within 0.01 r-soo- However, these fits 
are unweighted, so may be adversely affected by outlying points as 
their measurement errors are not accounted for. Instead, we directly 
fit the pedestal model of Equation |4] to determine if there is a sta- 
tistical reason to accept a pedestal addition to the profile. We find 
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Figure 7. Locally weighted regression fits to the scaled entropy profiles of 
the CC (white dotted line) and the NCC (black dashed line) groups. The 
confidence regions show the standard err ors on the fit. There is no evidence 
for the need of an entropy pedestal (e.g. iDonahue et al.ll20o3) to describe 
the CC profile at small radii. The top panel shows the number of groups 
contributing to the entropy profile fit as a function of radius for the CC 
(black) and NCC (grey) samples. The behaviour outside O.Srsoo is driven 
by a minority of systems. 



the following best-fitting relation when applying a weighted (with 
weights equal to the inverse measurement variance at each point) 
non linear regression fit, where the errors quoted are 1 a: 

/ T \ 0.79 (±0.04) 

S''(r) = 0.3 (±1.3) + 285 (±17) , (5) 

V ^500 / 

suggesting the pedestal modification is not statistically required, as 
5*0 is consistent with zero at the level of 1 a. Numerical simulations 
that employ only radiative cooling and the effects of gravity have 
shown the value of the entropy pedestal in CC clu sters to tend to 
zero as the system ages (Ettori & Brighe iitil 20081). In th e coolest 
system considered (4 keV) by Ettori & Brighentil (I2008"). this oc- 
curs in approximately 2 Gyr, and the entropy profile of the system 
is also seen to steepen during this time. Unfortunately, we do not 
have enough radial coverage for all groups to test the success of a 
pedestal model on individual group profiles. 

We can also see a suggestion of flattening in the entropy pro- 
files of CC groups at larger radii, between 0.5 rsoo and rsoo (see 
Figure |7](. There is also an apparent up-turn in the NCC profile at 
large radius. The top panel of Figure|7]shows the number of groups 
contributing at each radius. The effect seen in the CC profile at 
large radius is driven solely by the group RGH 80 as the radial cov- 
erage for the majority of CC groups extends only as far as 0.5 rsoo 
(five groups each add one data point at radii larger than 0.5 rsoo). 
Similarly, the apparent up-tum in the NCC profile is driven by only 
6 data points contributing outside 0.35 rsoo, and hence this feature 
should be regarded with caution. Therefore, comparing the profiles 
in Figure |7] the greatest difference between the CC and NCC pro- 
files occurs within 0.15 rsoo, where the NCC systems show higher 
entropy levels than the CC systems. At approximately 0. 15 rsoo the 



Figure 8. Power-law fits to the entropy profiles of the groups in the sample, 
divided into the CC (white dotted line with black 95% confidence region) 
and NCC (black dashed l ine with grey 95% confi dence region) groups. The 
theoretical slope of r^'^ iTozzi & Normanl200lb is also shown. 



profiles appear to converge, ignoring the artifacts introduced by the 
fits at large radius, as described above. 

As we are interested in comparing the slopes of the CC and 
NCC samples directly, we proceed in Figure [S] to fit power-law 
models in log-log space to the CC and NCC samples separately. 
These take the form. 



S'(r) = Sgn. ( — 



(6) 



V rsoo - 

which allows the derivation of the slope of the entropy profiles, 
allowing a comparison to that expected from gravitational effects 
alone. As already indicated in the locally weighted regression fits 
used in Figure [T] we see a significant difference in the observed 
slopes of the entropy profiles of the two samples. Figure [8] shows 
the 95% confidence bounds on the power-law fits, and at this level, 
the power-law fit to the NCC groups is flatter than that of the 
CC groups, within a radius of ~0.15 rsoo. The fitted slope to the 
CC systems is 0.71 ±0.02, compared to 0.57±0.04 for the NCC 
group s, both flatter than the expectation for gravitational effects 
alone jTozzi & Normaij200ll) . Outside ~0.15r5oo, the two power- 
law fits begin to converge, suggesting that the processes affecting 
the core have little effect at larger radii. As we will further show in 
Section |4.1.4| the NCC groups also show considerably more scatter 
in their entropies at small radii, evidenced by the broad confidence 
region in Figure|7] compared to their CC counterparts. 



4.1.4 Individual groups 

To look at the properties of individual groups, and how they relate 
to the mean trends, we fitted the scaled entropy profile for each 
group with a power-law model. In addition to the index of the fit, 
we determined the value of the entropy at three characteristic radii: 
0.01 rsoo, 0. 1 rsoo, and 0.5 rsoo- The values of each of these param- 
eters is shown in Table [3] We also fitted a power-law to the com- 



© 2009 RAS, MNRAS 0Q0.[Tll24l 



12 R. Johnson, T. J. Ponman and A. Finoguenov 



Table 3. The results of the power-law fits to the scaled entropy profiles 
of individual groups. Shown are the fitted index, and the scaled entropy at 
0.01 rsoo, 0.1 rsoo and 0.5 rsoo. The result for fitting a power-law to all the 
groups simultaneously is also shown. 
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bined data from all groups, which yielded a slope of 0.70±0.02. 
The slope and entropy measurements for each group indicate how 
each group varies from the mean, or expected theoretical, trends. 
NCC clusters are s een to show mor e scatter in their entropy pro- 
files at small radii jMcCarthv et~Z|[2008.') . This trend is also seen 
in our group sample, as for CC groups the rms scatter in scaled 
entropy at 0.01 rsoo is 6 keV^''^ cm^, whereas for NCC groups it 
is three times larger: 18 keV^'''' cm^. In contrast, at 0.1 rsoo, the 
profiles are much more comparable. The mean scaled entropy at 
0.1 rsoo for the CC groups is 70 keV^'''^ cm^, with an rms scat- 
ter of 26 keV^^"^ cm^, while for NCC systems the corresponding 
figures are 78 keV^''"^ cm^ and 23 keV^'^^ cm^. This lends more 
support to the observation in Section 14.1.3 1 that the CC and NCC 
entropy profiles converge outside 0.15 rsoo- 

The fitted parameters for Hickson 15 show it to have a rela- 
tively steep entropy profile, with very low central entropy. These 
are unexpected properties for NCC groups. Conversely, the fit to 
NGC 5 171 shows a very flat slope and high central entropy. How- 
ever, both these groups are current mergers, and hence likely to 
be far from equilibrium. N GC 5 171 has been studied in detail by 
lOsmond &. Ponmanl i2004) , and HCG 15 shows signs of having un- 
dergone a recent merger (E. O' Sullivan, private communication). 
Moreover in both these groups, the two-dimensional data cover 
only a small radial range, and the entropy fit has to be extrapolated 
well beyond the data to provide the values at the inner and outer 
radii tabulated in Table[3] Anomalous results for these two systems 
are therefore not unexpected. 
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Figure 9. A comparison of the pressure profiles between Cool (dotted white 
fine) and Warm (dashed black line) groups. The curves are locally weighted 
regression fits to the Cool groups and Warm groups separately. The confi- 
dence regions show the standard errors associated with the fit. The grey data 
points show individual scaled pressure measurements for the Cool groups 
(open circles) and the Warm groups (open squares). 



4.2 Pressure profiles 

Applying the modified entropy scaling of lPonman et alj(l2003h . the 
gas density scales as p oc T^^^, resulting in a scaling for pressure 
ofP (X T^/^To compare across the group sample, we therefore 
scaled the pressure profiles by T 

In an analogous manner to the analysis of the entropy pro- 
files, we again divided the groups into the Cool and Warm sam- 
ples, to see if there were any systematic differences in the pres- 
sure profiles. F-l-06 found that fitting a power-law to the pressure 
profiles di d not yield a good fit, and they instead followed the ap- 
proach of ISanderson, Finoguenov & Mohj | |2005 ) in fitting a lo- 
cally weighted regression curve to the data. This has the advantage 
of yielding a fit to the data that is independent of any model as- 
sumptions, and which also suppresses scatter. We follow the same 
approach here, and fit a non-parametric curve to the data (in log-log 
space) using the R function 'LOESS'. 

Figure |9] shows the fitted pressure profiles to the Cool and 
Warm group samples. The cooler groups are seen to exhibit steeper 
pressure profiles in the inner regions, most prominently within ap- 
proximately 0.05 rsoo, rising to a larger value of scaled pressure at 
the innermost radius. Outside ~ 0.15-0.2 rsoo, the fits to the two 
sub-samples begin to converge. It is possible that we are seeing 
the effect of the brightest group galaxy on the pressure profiles 
here, as the difference is most apparent within a small radius. The 
characteristic size of the BGG can be estimated using D25, the di- 
ameter of the isophote where surface brightness in the B-band is 
25 mag/arc sec^. We converted va lues of D25 from the RC3 galaxy 
catalogue dde Vaucouleurs et aljl 991) into physical radii using the 
distance information given in Tables[T]and|2] We find the mean size 
of the BGGs to be 0.05±0.02 rsoo- Separating this into the Cool 
and Warm samples yields the same result, i.e. there is no system- 
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Figure 10. A comparison of the pressure profiles between CC (dotted white 
line) and NCC (dashed black line) groups. The curves are locally weighted 
regression fits to the CC and NCC groups separately. The confidence re- 
gions show the standard errors associated with the fit. 



Figure 11. A comparison of the gas density profiles of Cool (white dotted 
line) and Warm (black dashed line) groups. The curves are locally weighted 
regression fits to the data, and the confidence regions show the standard 
errors on the fit. 



atic difference in the sizes of the central galaxies as a fraction of 
'"500- We postulate that the difference in pressure therefore arises 
from a difference in the dominance of the stellar mass component 
in the BGGs of Cool and Warm groups. 

We can examine this idea by comparing the masses of the 
group cores, and comparing with the luminosities of their central 
galaxies. Assuming hydrostatic equilibrium, the total group mass 
within radius r is given by. 



M(r) 



(7) 



Gp dr ' 

where P is the pressure and p is the density of the gas. We eval- 
uated Equation |7] for the Cool and Warm samples separately, us- 
ing the 'canonical' fits to the pressure (Figure |9} and gas den- 
sity profiles (see Figure [TTt. and removing the applied temperature 
scaling. Within 0.05 rsoo (roughly the optical extent of a typical 
BGG), the masses of the Cool and Warm systems are 2. 1 x Mq 
and 3.7xlO^^M0 respectively. Taking into account the differ- 
ence in rsoo between the Cool and Warm groups, these recovered 
masses are typical of brightest group galaxies at this radius (e.g. 
iHumphrev etai]|2006l) . 

We can test the stellar dominance by calculating the mass- 
to-light ratio of the Cool and Warm systems. The mean B-band 
BGG luminosities (with standard errors) for the Cool and Warm 
samples are 4.9±0.9x lO"^ Lb,o and5.7±0.9x lO"^ Lb,q respec- 
tively. Using the above masses then gives B-band mass-to-light ra- 
tios, within 0.05 rsoo, of 4.3 and 6.5 for the Cool and Warm groups 
respectively. This indicates increased dominance from the stellar 
mass in the cooler groups, which would lead to a more concentrated 
mass profile, and hence an excess pressure relative to the warmer 
groups. 

Dividing the sample into the CC and NCC classes, the same 
fitting procedure can be applied, and the fitted pressure profiles of 
the sample split on the basis of their core properties are shown in 



Figure [Tol The CC groups exhibit somewhat higher pressure than 
the NCC groups, although the two profiles overlap at ~ 0.15 rsoo- 
The typical radius of maximum temperature was found to be 
~ 0. 14 rsoo (see Section 13. 2t . so the excess pressure is seen only 
within the region of the CC. This is to be expected, since hydro- 
static equilibrium (equation 7) gives a steeper pressure gradient (for 
the same M{r)) when gas density is higher, as it is within CCs. 



4.3 Gas density profiles 

We also attempt to draw conclusions on the systematic differences 
between groups by examining the gas density profiles of the sam- 
ple. The modified entropy and pressure scaling applied to the group 
sample so far implies a scaling for gas density of p oc T^'^^, which 
we apply to the sample. 

We follow an analogous procedure to examining the pressure 
profiles of the systems in examining the gas density profiles of the 
Cool and Warm systems. We performed locally weighted regres- 
sion fits on the gas density profiles of the Cool and Warm groups 
separately, in log-log space, using the R 'LOESS' algorithm. The 
results of these fits are shown in Figure[TT] The gas density profiles 
of the Warm and Cool groups again appear to diverge at approx- 
imately 0.05 rsoo, with the Cool groups showing a higher scaled 
central density compared to the Warm groups. As explained in 
Section l42l this is showing the effect of the central galaxy. The 
pressure of the Cool systems is more strongly peaked, and conse- 
quently, the gas density is affected in the same manner. 

The same comparison can be made dividing the sample on the 
core properties of the groups, and this is shown in Figure[T2] In line 
with the trends seen in the pressure profiles, the scaled density pro- 
files of CC groups are higher than those of NCC groups, within a 
radius of ~ 0.15 rsoo - This is showing the effects of radiative cool- 
ing in the CC systems, leading to a higher density and pressure in 
the region of the core (see Section l4~2] l. 
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Figure 12. A compaiison of the gas density profiles of CC (wfiite dotted 
line) and NCC (black dashed Hne) groups. The curves are locally weighted 
regression fits to the data, and the confidence regions show the standard 
errors on the fit. 



5 ENTROPY DISTRIBUTIONS 

Considerable insight into the heating and cooling processes affect- 
ing the intracluster medium can be gained from looking at how the 
entropy is distributed as a function of gas mass of the system (e.g. 
IVoit et al" 2003). Having fitted power-laws to the scaled entropy 
profiles of all the groups individually (see Section |4.1.4t . we used 
these alongside /3-model fits to the gas density profiles to determine 
the gas mass in a series of radial shells. The shells were stepped in 
radius from 0.0 to 0.5 ^500, with width 0.0005 r5oo, and we set the 
characteristic radius for a shell to be its mean radius. The fits to the 
entropy and gas density profiles were used to interpolate the en- 
tropy and density at the these characteristic radii. The gas mass in 



each shell, Mgas was determined as follows, 



■ A-K r p Sr, 



(8) 



where r is the characteristic radius of the shell, p is the density eval- 
uated at radius r, and 5r is the width of each shell, i.e. 0.0005 rsoo. 

With entropy and gas mass then described by the same set 
of radial bins, we set up a series of scaled entropy bins of width 
15 keV^''^ cm^, and determined the total gas mass in each entropy 
bin. To be able to compare groups directly, we scaled the gas mass 
by P to remove the dependence of gas mass on the mass of the 
system as a whole. This is a result of the applied entropy scaling 
which gives a gas density scaling oc T^^^ (see Section \43\ . In a 
particular shell of mass Mgas, 



— A-K r 



3 

r p. 



From the virial theorem, the radius r is proportional to T 
with the density p scaling as T ^^'^ we find. 



(9) 
and 

(10) 



Dividing the sample into CC and NCC systems, we can ex- 
amine the entropy distribution in both types of system. Figure [T3] 
shows a histogram of the gas mass as a function of entropy. This 



Figure 13. Histogram of gas mass as a function of entropy for the CC 
groups, with each group colour-coded as shown in the legend. This shows 
the total gas mass in a particular entropy bin, indicating where the majority 
of the gas mass lies. We have used the modified entropy scaling of T^/^ to 
directly compare the systems, and have scaled the gas mass by (see text 
for details). 
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Figure 14. Histogram of gas mass as a function of entropy for the NCC 
groups, with each group colour-coded as shown in the legend. This shows 
the total gas mass in a particular entropy bin, indicating where the majority 
of the gas mass lies. We have used the modified entropy scaling of T'^/^ to 
directly compare the systems, and have scaled the gas mass by f ^ (see text 
for details). 
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indicates where the majority of the gas lies in terms of its entropy 
value, for each of the CC groups. The shape of these histograms re- 
quires careful interpretation. Since entropy rises outwards, the peak 
in these histograms is produced as a consequence of the limiting 
upper radius of the gas mass calculation (0.5 rsoo). At this radius, 
some groups have reached a higher entropy level than others, which 
produces the high entropy tail in the plot. It is evident is that there 
are five CC systems (HCG 42, NGC 2300, NGC 2563, NGC 5846 
and SRGbl 19) where higher entropy values are reached at 0.5 rsoo 
compared to the remainder of the CC sample. NGC 2300 reaches 
very high entropies (> 600 keV^''^ cm^), although the entropy pro- 
file is extrapolated from just 4 data points, all within 0.2 rsoo, so 
the result shown here should be interpreted with caution. The CC 
systems generally show a sharp increase in scaled gas mass with 
scaled entropy on the leading edge of the peak, and show a dom- 
inant peak at approximately 105-120keV^'^'^cm^, indicating that 
similar entropy levels are reached by the groups at 0.5rsoo. 

The entropy distributions for the NCC systems are shown in 
Figure[T4] In comparison to the CC systems, the NCC groups show 
a wider diversity in entropy distributions, as can also be seen in 
Figures [15] and [16] This is in agreement with the observed tight 
trends for the entropy profiles of CC groups, seen in Figure[8] whilst 
the NCC systems show a wider spread in their entropy properties. 

The entropy histograms also show a lack of lower entropy gas 
in the NCC systems — in the CC systems the mean scaled gas mass 
with scaled entropy < 90 keV^^^ cm^ is ~9.8xl0" Mq keV^, 
compared to ~5.OxlO"'M0 keV~^ in the NCC systems, thus 
confirming the results seen in Section 14.1.31 Also, it is apparent 
that the entropy at which scaled gas mass peaks is slightly different, 
with the CC systems peaking in the bin 105-120 keV cm^, whereas 
the NCC histogram peaks at 120-135 keV cm^. 

Since this analysis employs scaled entropy and scaled gas 
mass, the mean effect of system mass should have been scaled out. 
We find the mean total scaled gas mass, evaluated within 0.5 rsoo, 
for the CC and NCC groups to be (4.7±0.5)x 10"Mq keV^ 
and (4.1±0.5)x lO"Af0 keV~^ respectively. Therefore, there is 
no significant difference in the overall scaled gas content of these 
systems within 0.5 rsoo, indicating that the different core properties 
do not result from global differences in gas content. 



5.1 Comparison to a theoretical entropy distribution 

In an effort to ascertain the magnitude of the effects of feedback 
processes acting on the observed entropy distributions, we can 
compare these distributions to theoretical distributions where the 
effects of non-gr avitational proces ses have been ignored. The ana- 
lytical models of lVoitetal.ll2003b show that the radial distribution 
of entropy is approximately proportional to the enclosed gas mass if 
the gas accreted is cold and t he acc retion is smooth. Looking specif- 
ically at groups, IVoit et al.l Jioosb postulate smoother accretion in 
groups relative to clusters, to explain the observed differences in 
the entropy profiles of systems spanning this range in mass. We can 
begin to evaluate the effects of feedback on the group gas by first 
computing the expected entropy histogram for a system where only 
gravitational processes are considered. The observed difference be- 
tween this 'theoretical' distribution and the observed distribution 
can then be used as a probe of the impact of feedback processes. 

We firstly require a baseline entropy profile S{r) for a system 
unaffected by non-gravita tional processes. For this , we adopt the 
baseline entropy profile of lVoit, Kay & Brvanl ( l200^ , derived from 
smoothed particle hydrodynamics (SPH) cosmological simulations 



of galaxy clusters, without any cooling or feedback processes. This 
has the form 

^(r-) = 1.32 S200 f — V'\ (11) 
V r200 / 

where 5200 is defined as 



5*200 ~ 362 keV cm' 



/T200 \ 



(12) 



\Tx J 

where we have dropped the cosmological factors from IVoit et al.l 
(2005), which are negligible given the low redshift of our sam- 
ple. For consistency within this work, we us e S to denote entropy , 
equivalent to the K used bv lVoit et al] ( |2002|) and lVoit et 



te e ntropy , 
ai]( l2005l) . 



For the temperature Tx we use our mean X-ray deriv ed temper- 
atures , and calculate r2oo in the manner described by IVoit et all 
( l2005h . This requires knowledge of both r2oo and the mass en- 
closed at this radius, M200, which is simply the volume at r2oo mul- 
tiplied by 200 peril - The former can be calculated by determining a 
con version factor from the known r spQ. We assume an NFW profile 
Navarro, Frenk & Whitell997h . evaluated at an overdensity of 



Sun et al. 1 2008h found a concentration at this overdensity of 



(e.g 
500. 

4.2 for their sample of galaxy groups observed with Chandra and 
we proceed to use this concentration in the following analysis. We 
find for an NFW dark matter halo with this concentration, that rsoo 
= 0.66 r2oo. This modifies Equation [TT] to the following, 

'0.66r^ 



S{r) = 1.32 520 



rsoo 



(13) 



Therefore, we use Equation [13] to describe the theoretical entropy 
profiles of our systems. 

T o calculate the gas mass, we assume, in line with lVoit et al.l 
( |2002[) , that the gas density profile follows the dark matter den- 
sity profile, s caling the latter by a factor Qt I f^dm, assuming 
= 0.022 h"^ (' Voit et afllioosl) and fldm = 0.2. We then apply the 
same analysis procedure as for the observed entropy distributions, 
calculating first gas density and entropy as a function of radius, be- 
fore determining the gas mass in each bin of entropy. We use scaled 
entropy and gas mass as for the observed histograms throughout our 
comparisons. Figures[T5]and[T6]show the observed and theoretical 
entropy distributions for the CC and NCC groups respectively. The 
rightmost bin of each distribution should be treated with caution, 
since due to our imposed cut-off at 0.5 rsoo, these bins are incom- 
plete. 

In all cases, the observed maximum entropy far exceeds 
the theoretically derived maximum entropy, suggesting substantial 
modification must have occurred to the theoretical distribution to 
yield the observed distribution. We now explore some simple heat- 
ing and cooling prescriptions, to see to what extent they are able to 
reproduce the modified entropy distributions observed. 

5.2 Entropy modifications 

IVoit et all ( |2002|) consider three simple modifications to a baseline 
entropy distribution, which we summarise below: 

(i) A truncation in the entropy distribution where gas is removed 
if its entropy falls below a threshold value corresponding to a given 
cooling time. This approximates the elfects of cooling, where gas 
cools out to form stars, or is heated by feedback from proximate 
cooled gas (e.g. via supernova feedback), raising its entropy such 
that it convects to larger radii. 

(ii) A shift in the entropy distribution, such that the entropy of all 
gas is increased by a certain baseline entropy, mimicing the effects 
of pre-heating. 
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Figure 15. The observed entropy distributions of the CC groups (shown in grey) and the calculated theoretical entropy distribution (shown in black). In all 
cases the observed distribution extends to a higher maximum entropy, compared to the theoretical distribution. The dashed line shows the entropy distribution 
resulting from a modification in the form of an entropy shift, followed by radiative cooling (see Section l53t . Note the extended x axis in the case of NGC 2300, 
to show the high entropies reached by this group (this is due to a significant extrapolation, see Section |5). The groups are ordered in increasing mean 
temperature, from top left to bottom right, and all curves extend to 0.5 rsoo. 



(iii) Lowering of the whole entropy distribution due to radia- 
tive cooling. The entropy distribution of the gas to the 3/2 power 
(i.e. S^^^) is reduced by the 3/2 power of the critical entropy (Sc) 
across the whole entropy distribution. Gas which drops below zero 
entropy as a result of this modification is removed, as in model (i). 
The form of the entropy reduction is based on the approximation 
that the cooling function for gas in the group regime (T < 2 keV) 
has a temperature dependence A(r) oc T~^^'^ . 

The critical entropy Sc. in models (i) and (iii) is some fraction of 
'S'2 00, the e ntropy at r2oo • In fact, following the discussion presented 
bv lVoit et a l. (2002), we specify, 



calculating Sc for each group, where T is its mean temperature. 

It can immediately be seen by considering the differences be- 
tween the observed (grey) and theoretical (solid black) distribu- 
tions in Figures [T5] and [16] that simple truncation (model (i)) and 
shift (model (ii)) models are inadequate. Truncating the theoretical 
distribution at a critical entropy would lead to a loss of low en- 
tropy gas, but would not modify the shape of the distribution above 
the critical entropy, making no improvement in the agreement be- 
tween the observed and theoretical distributions at high entropy, 
and making things worse at low entropy (where low S gas is ac- 
tually observed in nearly all groups). If instead we add a constant 
(a 'shift') to the theoretical entropy distribution, we simply move 
this distribution in comparison to the observed distribution. This 
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Figure 16. The observed entropy distributions of the NCC groups (shown in grey) and the calculated theoretical entropy distribution (shown in black). Like 
the CC systems, in all cases the observed distribution extends to a higher maximum entropy compared to the theoretical distribution. The dashed line shows 
the entropy distribution resulting from a modification in the form of an entropy shift, followed by radiative cooling (see Section|5]2). The groups are ordered 
in increasing mean temperature, from top left to bottom right, and all curves extend to 0.5 rsoo . 



will improve the agreement at high entropy, but would leave the 
distribution devoid of any low entropy gas. 

The radiative cooling modification implemented bv lVoit et al.l 
( I2OO2I) reduces the theoretical entropy distribution S^^^ by Sj 
and removes all gas with an entropy less than zero. The impli- 
cations of such a model are clear; unlike the truncation and shift 
models, low entropy gas will remain. However, it is not possible 
to increase the maximum entropy reached by the theoretical model 
in this example, as whatever the value of Sc, the resulting entropy 
is a reduction of the original. Therefore, a simple radiative cooling 
model does not aid in matching the entropy profiles at the highest 
entropies reached in Figure s[T5] and [T6] 

It seems clear that a more sophisticated approach is required, 
and we can envisage a model which combines the simple modi- 
fications described above. To fix the high entropy behaviour, we 
first need a shift in entropy, which we allow to subsequently cool 
via the radiative modification to populate the low entropy end of 
the distribution. In physical terms, this model uses an entropy shift 
to mimic early pre-heating of the group gas, which subsequently 
cools through radiative cooling. In fact, we can calculate the ex- 
act level of the required entropy shift, S shift, for the final maxi- 
mum theoretical entropy to match the maximum entropy achieved 
by the observed distribution. Assuming the maximum observed en- 
tropy So,max is achieved by first adding a shift Ssm ft to the max- 
imum theoretical entropy St, max, before subsequently cooling as 
described above by an amount Sc^'^, we have, 

So!max ~ {St, max + Sshift)^ ^ ^ ~ Sj . (15) 



Rearranging Equation [15] we have 

Sshift = {So^max + 'S'^'' — St,max- (16) 

Therefore, applying the shift Sghift before cooling by Sc^^ will 
match the right hand edge of the entropy distribution. The dashed 
lines in FiguresfTSlandl I6lshow this entropy modification to the the- 
oretical distribution. In the case of NGC 2300, we have extended 
the X axis to show the high entropy behaviour at 0.5 rsoo (see Ta- 
ble |4X4ll, a result of fitting a power-law to a profile of only four 
data points, all located within 0.2 rsoo - Therefore, due to the exten- 
sive extrapolation the entropy distribution of NGC 2300 should be 
treated with caution. 

Our applied "shift -I- cool" model is certainly more successful 
than models (i)-(iii), but it is able to produce a reasonable represen- 
tation of the observed distribution in only a small number of cases. 
In the CC sample, NGC 4636 and RGH 80 are best represented by 
the applied model, however, Hickson 97 and SS2bl53 also perform 
reasonably well. In the NCC sample, Pavo, Abell 194, NGC 5I7I 
and 3C 449 all show a reasonable agreement between the observa- 
tion and the applied model. However, the vast majority of groups 
show marked differences between the "shift -I- cool" model and the 
observed entropy distribution. In every case this disagreement is 
due to a lack of low entropy gas in the model. Conversely a higher 
gas mass than observed is usually present in the highest entropy 
bins. One caveat with our model is that the total gas mass within 
0.5 rsoo remains fixed. In reality, this value would change, due to 
movement of the gas as a consequence of the applied heating and 
cooling mechanisms. For example, if low entropy gas is removed 
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from the distribution, higher entropy gas will flow in. Boosting the 
entropy of the gas would in reality cause it to expand, such that 
a proportion of the gas mass at high entropy would actually fall 
outside O.Srsoo, and the applied entropy shift will have been un- 
derestimated. Therefore, once the expansion of the gas is taken into 
account, the applied entropy shift needs to be somewhat larger than 
shown in Figures[T5land[76l This will exacerbate the difficulties of 
the model, which arise (see discussion below) from the large en- 
tropy boost required in many systems. 

The required entropy shifts range between llSkeVcm^ 
(RGH 80) and 567keVcm^ (NGC 2300). The models which re- 
create the observed entropy distribution most successfully are those 
where the value of the applied entropy shift is very similar to the 
critical entropy Sc (Equationll4t which governs the radiative cool- 
ing (Equation |15t , since in these cases, the shift provides a good 
match at the high entropy end of the distribution, whilst the cool- 
ing is able to repopulate the low entropy end. For example, in 
NGC 4636, the applied entropy shift is 148keVcm'^ compared to 
a critical entropy of 164 keV cm^ and in RGH 80, the entropy shift 
is 1 15 keV cm^ compared to a critical entropy of 109 keV cm^ . The 
groups for which this works are those with the smallest difference 
in maximum entropy between the observed and theoretical distribu- 
tions. In the majority of cases however, the required shift to match 
the high entropy end is much larger than the critical entropy, in 
which case the gas entropy is boosted to the point where cooling 
has little effect. 
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Figure 17. The required entropy shift S^y^i fi versus the critical entropy Sc 
for the CC (open circles) and NCC (filled squares) groups in the sample. 
The solid fine is the line of equality, and the dashed fine shows the same but 
re-normahsed by the factor '~1.6. 



5.2.1 The effects of metallicity 

The lVoit et al.l ( |2002|) prescription for the critical entropy assumes 
that cooling has proceeded for 15 Gyr, and that the emissivity of 
the gas corresponds to a metallicity of 1/3 solar. However, in the 
group regime, emissivity is a strong function of metallicity (due 
to strong line cooling) and in CC groups there is a strong central 
abundance gradient (Rasmussen & Ponman 2007, Johnson et al., in 
prep.) leading to central metallicity close to solar values. We can 
assess the potential effects of such an abundance gradient on our 
'shift -I- cool' model by instead assuming solar metallicity. This is 
assumed across the radial range, so will in fact be an over-estimate 
of the effects of a metallicity gradient. For a hot plasma at 1 keV, 
the emissivity scales up by a factor of ~ 2.0 wh en the metallicity 
increases from 1/3 solar to solar. Equation 17 of IVoit et al.l j2002h 
describes the rate of change of entropy as follows, 

^ oc r^''^A(r), (17) 

where our S is equivalent to the lVoit et alj ( |2002|) K. The lVoit et al.l 
( l2002h model assumes the right hand side of this relation to remain 
constant, such that it can be simply integrated, 

j dS^'"^ Qc J T^^^A{T)dt. (18) 

Integrating over a Hubble time (tn) yields, 

Sc oc A(r)fH]'^', (19) 

Hence the effect of a change in emissivity, caused by an increase 
in metallicity, is to increase the value of Sc. The factor increase 
is equal to the factor increase in emissivity, to the power 2/3 (i.e. 
2.0^^/^^ « 1.6). InFigureini 

we show the size of the required en- 
tropy shift Sshi ft versus the critical entropy Sc in our 'shift -I- cool' 
model. Since almost all groups contain observable low entropy gas, 
the model can only be successful if Sc is similar to or larger than 



Sshift- For the unmodified cooling threshold (solid line) it can be 
seen that only a few groups satisfy this requirement. However, with 
the factor 1.6 scaling in Sc, the model could be successful for at 
least half the groups, although approximately one third of the sys- 
tems still require entropy shifts much greater than the critical en- 
tropy. 

Our model is still very simple, of course. The uniform solar 
abundance will tend to overestimate the effects of cooling, as will 
the generous timescale of 15 Gyr allowed for the cooling to pro- 
ceed. On the other hand, we have not allowed for the hierarchi- 
cal merger history of the systems. This requires investigation with 
more sophisticated cosmological models. Nonetheless, we tenta- 
tively conclude that the magnitude of the entropy shift required in a 
significant fraction of groups is so high that our 'shift + cool' model 
will be unable to reproduce their properties. 

5.3 Effects of feedback 

We turn now to the link between modification of the entropy dis- 
tribution and the feedback processes which are most likely re- 
sponsible for this modification. If we make the assumption that 
the most prominent source of feedback in galaxy groups comes 
from the member galaxies, either through energy input from su- 
pernovae or feedback from active galactic nuclei (AGN), we can 
begin to quantify the likely impact of feedback. We construct a 
rough measure of 'feedback impact' using the ratio of the op- 
tical luminosity of the group to the thermal energy of the gas: 
fgrp = L B. group /{MgaaTx), where Is the mean X-ray tem- 
perature of the system, and Mgas is the total gas mass at 0.5 rsoo 
(derived as described in Section[5). The motivation here is that inte- 
grated feedback from supernovae should sc ale approximately with 
the stellar mass, and given the relationship jMagorrian et al.|[l998l) 
between the mass of supermassive black holes and the stellar mass 
of the spheroids they inhabit, integrated AGN feedback might also 
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Figure 18. The entropy distribution for stacked CC groups (top row) and stacked NCC groups (bottom row), in three bins of feedback impact, fgrp (see text 
for detinition). Feedback impact is increasing from left to right in both rows, from 0.0 < fgrp ^ 0.22 to 0.22 < fgrp ^ 0.6 and the rightmost panels show 
fgrp > 0.6. As the impact of feedback increases, so to does the maximum entropy reached by the groups, both for CC and NCC systems. 



be expected to scale roughly with stellar mass. The effect of a given 
amount of injected energy on the intergalactic gas will depend on 
the factor by which it increases the existing total thermal energy, 
which motivates the denominator, (MgasTx), in our expression 
for fgrp- This measure of feedback impact therefore has units of 
Lb,q Mq~^ keV~^. For clarity we will not show the units in the 
remainder of this section. 

Our measure of feedback impact is admittedly imperfect. In 
the first place, the A'-band luminosity would provide a better mea- 
sure of stellar mass, since the B-band is strongly affected by young 
stellar populations. A full stellar population analysis would be even 
better, of course. However, we have consistent measurements of 
B-band luminosities for all the groups in our sample, and are en- 
couraged to believe that biases due to young stars will not have a 
major impa ct, by the fact that X-ray b right groups have a low spi- 
ral fraction faelsdon & Ponmanll2003h . A second respect in which 
our formula for fgrp is imperfect is that the denominator actually 
reflects the thermal energy of the gas after the feedback has had its 
effects, rather than before. We explore the effects of modifying the 
denominator at the end of this section, and find that it makes no 
qualitative difference to our conclusions. 

Assuming that fgrp is monotonically related to feedback im- 
pact, we now use it to divide our sample into 3 subsets of (almost) 
equal size, covering ranges in feedback impact of 0.0 < fgrp ^ 
0.22, 0.22 < fgrp ^ 0.6 and fgrp > 0.6. These bins contain 10 (8 



CC + 2 NCC), 9 (6 CC -H 3 NCC) and 9 (4CC + 5 NCC) groups 
respectively. Figure [TS] shows the stacked entropy distributions for 
the CC (top row) and NCC (bottom row) groups, shown in order of 
increasing feedback impact from left to right. In both samples, we 
see a trend for an increase in the amount of gas at higher entropies 
as the feedback impact increases. There is also a suggestion of the 
peak of the distribution moving to higher entropies as the level of 
feedback impact increases. Qualitatively, this appears to indicate 
that the gas in groups with the largest capacity for feedback from 
their galaxy members (in the form of energy deposited into the in- 
tracluster medium from supernovae explosions or AGN) has been 
pushed to higher entropies. This is true regardless of whether the 
system has a CC or not. 

Given the similarity of the entropy distributions in the highest 
feedback bins, we choose to combine these bins. Therefore, the 
sample now consists of a 'lower' feedback bin (fgrp ^ 0.22) and a 
'higher' feedback bin (fgrp > 0.22). Table|4]shows the break-down 
of our groups into the CC/NCC and high and low feedback classes. 
In brackets we show the expected number of groups in each of the 
four categories, given the relative ratios of CC to NCC groups, and 
high to low feedback systems (18 to 10, in each case) if the two 
classifications (CC/NCC and high/low feedback) were completely 
independent. This suggests that there is some tendency for the NCC 
systems to have a higher feedback impact than CC systems. The 
median feedback impact for the whole sample is 0.31, whereas for 
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Table 4. The distribution of CC and NCC systems between the higher and 
lower feedback bins. The numbers in brackets show the expected number of 
groups in each of the four categories, if the classifications of CC/NCC and 
high/low feedback are independent. 







CC 


NCC 


Total 


High Feedback {fgrp 


^ 0.22) 


10(11.6) 


8 (6.4) 


18 


Low Feedback (fgrp 


< 0.22) 


8 (6.4) 


2 (3.6) 


10 


Total 




18 


10 





CC systems it is 0.27 and for NCC systems 0.62. We quantify the 
significance of the relationship between feedback and CC status by 
calculating how likely it is that there should be at least 8 groups 
in both the high feedback, NCC class, and the low feedback, CC 
class, given that only 6.4 systems were expected in each, using the 
binomial probability distribution. This probability is 6 percent, so 
we conclude that there is some evidence for a relationship between 
the two properties, but it is clear from Table |4]that the relationship 
is rather weak. 

In addition to affecting the energy of the gas in the intracluster 
medium, feedback processes can also inject and redistribute met- 
als in the intracluster medium. Supernova feedback will introduce 
more metals into the intracluster medium, whilst AGN will not, so 
by examining the metal content of the gas in high and low feedback 
systems, we may hope to differentiate between the two sources of 
feedback. Figure[79lshows the metallicity profile of the lowest feed- 
back {fgrp < 0.22) systems and the higher feedback systems (fgrp 
^ 0.22). We have stacked the abundance profiles of the two feed- 
back sub-samples into equal radial bins. The lowest feedback sys- 
tems extend to a smaller radius than the higher feedback systems. 
We therefore show the contribution at small radii from the lowest 
feedback systems as a single grey point. Over the majority of the 
profile, the low and high feedback systems are remarkably similar. 
Only in the innermost radial bin do the profiles differ significantly, 
with the higher feedback systems showing lower central metallic- 
ities than the lowest feedback systems. We estimated an approxi- 
mate integrated metal fraction for each group by determining the 
product of metallicity and gas mass, summed over a series of radial 
shells, and dividing by the total gas mass. We restricted this calcu- 
lation to within 0.3 r^oo, to ensure that the majority of groups have 
the appropriate radial spectral coverage to avoid significant extrap- 
olation of the metallicity profile. For the low feedback systems, we 
find the mean metal fraction (with standard error) to be 0.23±0.04, 
whereas the mean metal fraction (with standard error) of the high 
feedback systems is 0.25±0.03. Hence, across this radial range, the 
integrated metal fraction of the two subsets is essentially identical. 
This is the case despite the abundance difference in the innermost 
radial bin since the metal fraction is dominated by shells at larger 
radii, where the majority of the gas lies. These results run counter 
to naive expectations. We would expect higher feedback systems to 
have experienced more supernova feedback, and therefore to show 
higher abundances and integrated metal mass fractions. 

The fact that we see no increase in the integrated metal mass 
fraction in the highest feedback systems, coupled with the lower 
central metallicity of these groups, suggests that AGN may provide 
the dominant source of feedback, as supernova feedback would in- 
ject metals into the gas, along with the deposited energy, leading 
to extra enrichment. We have already seen in Figure [TS] that there 
is a lack of low entropy gas in the high feedback systems, com- 
pared to the low feedback groups. We can examine whether high 
metallicity, low entropy gas has been boosted in entropy (generat- 
ing high entropy gas with high metal abundance) by considering 



the distribution of metals as a function of entropy, rather than ra- 
dius. In Figure[20l we show the stacked abundance as a function of 
entropy for the low {fgrp ^ 0.22) and high {fgrp > 0.22) feedback 
systems. The samples have been divided into the same bins in en- 
tropy for comparison. As the low feedback systems extend to lower 
entropy, we show the remaining data points outside the range of 
comparison as an open grey circle in Figure |20] In the highest en- 
tropy gas, the gas in the low and high feedback systems is similarly 
enriched. The lowest entropy gas in the highest feedback systems 
is of lower metallicity, compared to gas at the same entropy in the 
lowest feedback systems. Hence, if in the high feedback systems, 
the high metallicity, low entropy gas has been boosted in entropy, 
it would need to be pushed beyond 0.5 rsoo for us not to detect it in 
our combined data. 

One concern in studying the relationship between our feed- 
back impact parameter and gas entropy is that both are related to 
gas density. Systems with low gas density will have high entropy, 
and will also (via the entry of Mgas into our expression for fgrp) 
tend to have high feedback impact. As mentioned earlier, the de- 
nominator of the expression for fgrp should really be the total ther- 
mal energy in the gas before the action of feedback, rather than 
afterwards. Assuming that all systems would contain a cosmic ra- 
tio of gas to dark matter, in the absence of cooling and feedback, 
we can construct an alternative measure of feedback impact, as 
Lb. group/ {MtotTx). For self-similar systems at a given epoch, the 
total mass and characteristic temperature are related by Mtot oc 
y3/2 therefore explore a modified definition of feedback im- 
pact of LB.group/Tx^"^ ■ 

Using this revised parameter to partition our group sample into 
high and low feedback classes, we find results very similar to those 
obtained with our earlier definition of feedback impact: gas with 
entropy greater than ~30 keV'^^^' cm^ is similarly enriched in the 
low and high feedback systems, and in the innermost entropy bin, 
the highest feedback systems reach lower central metallicity. We 
also find a trend for the higher feedback systems to reach higher 
entropies than the lowest feedback groups, as seen in FigurefTSjfor 
our unmodified measure of feedback impact. We conclude that the 
gas density dependence of our original feedback measure does not 
affect the results. 



5.4 The o -T relation 

We can further investigate the impact of feedback by considering 
the quantity Pspec, which is the ratio of the specific energy of the 
galaxies to the specific energy of the gas, 

Psp.. = i^, (20) 

where fj, is the mean mass per particle, nip is the proton mass, k 
is Boltzmarm's constant, and T is the temperature of the gas. If 
feedback raises the energy of the gas, then the consequence is to 
reduce (3spec, as the subsequent specific energy of the gas would 
dominate over the specific energy of the galaxies. 

We can use the a — T relation to investigate where groups lie 
with respect to fUspec = 1- From a theoretical perspective we expect 
to see similar specific energy in gas and galaxies, unless feedback 
has had a major effect, and in the case of rich clusters, observations 
have shown that relaxed systems generally lie along a trend in the 
a — T plane of the form a ocT'^'^, whilst g roups appear to follow 
a steeper trend (see the discussion in lOsmo nd & Ponman (2004)). 
We have used a BCES orthogonal regression ( lAkritas & Bershadvl 
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Figure 19. The stacked abundance profiles of the groups with the lowest 
level of feedback impact (fgrp < 0.22, see Section [53), shown as black 
filled squares and the groups with higher levels of feedback (fgrp ^ 0.22), 
shown as red filled triangles. Each sample has been divided into equal bins 
of radius. Vertical error bars show the standard error in each abundance 
bin, and horizontal error bars show the width of each radial bin. The radial 
coverage of the low feedback systems extends to a smaller radius than the 
high feedback systems, and we show these remaining low feedback systems 
as a single grey open circle point. 



Figure 20. Stacked abundance for the low feedback (black filled squares) 
and the high feedback (red filled triangles) systems shown in Figure [T9l 
but now as a function of scaled entropy. Each sample has been divided into 
five equal bins in scaled entropy. Vertical error bars show the standard error 
in each abundance bin, and horizontal error bars show the width of each 
entropy bin. The lowest feedback systems extend to a lower entropy, and 
we show these remaining low feedback systems as a single grey open circle 
point. 



1 1996b to fit a straight line to the cr — T relation for our group sample, 
yielding the fit: 

logio a = 0.9 (±0.3) logw T + 2.57 (±0.02). (21) 

This method accounts for both the measurement errors and the un- 
known intrinsic scatter in the relation. 

Figure |2T] shows the a — T relation for the group sample, 
where we show the line corresponding to Pspec = 1 as a solid 
line. Six of the ten NCC systems are seen to lie above the line of 
Pspec ~ 1, whereas the majority of CC systems lie beneath this 
line. If the lack of central cooling in NCC system s resulted from 
strong feedback, as in the model o f lMcCarthv elaP |2008), who in- 
voke pre-heating prior to cluster collapse, then the specific energy 
of the gas should be higher than that of the galaxies in NCC groups, 
i.e. /3spec < 1, whereas in practice, we see the exact opposite. This 
may indicate the importance of cluster merging in the formation of 
NCC systems, as in merging systems, the velocity dispersion can 
be boosted, which would in turn increase Pspec- This is supported 
by the NCC group with the highest Pspec, Hickson 92 (Stephan's 
Quintet), wh ich is clearly a system currently undergoing multiple 
mergers (e.g. lTrinchieri et alj2003l) . 

5.5 Discussion 

We can connect the evidence of Section|5]to build a coherent picture 
relating galaxy group properties, entropy and the effects of feed- 
back. Section lS^ shows that modifications must have been made to 
the gas entropy distribution in these groups. Comparing with a the- 
oretical distribution, which does not allow for any non-gravitational 
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Figure 21. The cr-T relation for the groups in the sample. Black open circles 
denote the CC groups, and grey filled squares denote the NCC systems. The 
solid line shows the relation (3spec = 1 , and the dashed line shows a BCES 
orthogonal regression fit to the data (see text for details). 
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processes, we find a higher maximum entropy in the observed en- 
tropy distributions. Simple entropy modifications are not sufficient 
to convert the theoretical distribution to the observed distribution, 
and we rule out simple entropy shift, truncation and radiative cool- 
ing models. Our model which mimics pre-heating followed by a 
bout of radiative cooling (our so-called 'shift + cool' model) only 
performs well in cases where the required entropy shift is modest. 
In many cases, the required large entropy shift suppresses cooling, 
making it impossible to populate the low entropy end of the ob- 
served distribution. 

To avoid this problem, we require a model which ap- 
plies a larger entropy shift to higher entropy gas. Two physical 
mechanisms which might achieve th is are entropy amplification 
jPonman et Zll2003l; pN^it et al.ll2003h or episodic heating. In the 
former, when pre-heated gas hits an accretion shock there is an 
entropy boost, and higher entropy gas hitting the accretion shock 
gets a bigger boost than lower entropy gas. In the case of episodic 
heating, small bursts of heating could raise the group entropy, and 
subsequent cooling operates more efficiently in the lower entropy 
gas. Hence, over a number of cycles, the highest entropy gas could 
be boosted more than the lower entropy gas, as required to model 
our observed entropy distributions. The former mechanism could 
require up to two orders of magnitude less energy than the latter, 
depen ding on when the pre-heating takes place HMcCarthv et al.l 
l2008h . 

Defining feedback impact as Lb, group/ (MgasTx), we find 
gas of higher entropy in systems with higher feedback impact. 
However, there is no significant increase in metallicity in the high- 
est feedback impact systems, and the metal mass fraction in the gas 
is similar to the lowest feedback systems. As our measure of feed- 
back impact scales with the optical luminosity of the group, we 
would expect an increase in the number of supemovae, and hence 
an increase in the metallicity of the highest feedback systems. The 
lack of such a trend suggests that metals must have been removed 
from the gas within 0.5 rsoo either by ejecting it to larger radii, 
or by locking up a larger fraction in stars or stellar remnants, and 
strongly favours AGN, rather than supernovae explosions, as the 
dominant source of feedback. 

We see higher metallicity gas in the centres of the lowest feed- 
back groups, but across the majority of the radial range, the metal- 
licity profiles of the lo w and high feedba ck systems are compara- 
ble. The conclusion of Ijetha et alj ( |2007[) that central AGN within 
groups influence the gas properties of groups on a local, rather than 
a global, scale is supported by this observation, as the effect of 
feedback is notice able only at the innermost radius. This con trasts 
with the results of ICroston. Hardcastle & Birkinshawl feOOSh . who 
found evidence for radio sources affecting the global properties of 
the group gas, as represented by their location in the luminosity- 
temperature plane. 

In the highest entropy gas we study, whose entropy exceeds 
that expected due to gravitational collapse alone, similar levels of 
enrichment are reached in high and low feedback syatems. This 
again points to an AGN origin for the bulk of the feedback, though 
this feedback may well have taken place within precursor struc- 
tures, rather than being driven by a central AGN within the assem- 
bled group. 

Eight of the ten NCC systems have feedback impact >0.22, 
putting them in our highest feedback bin, and five of these have 
a value of (3spec > 1. This is puzzling, since high feedback im- 
pact should indicate raise the specific energy of the gas, resulting 
in Papec < 1. The observed high Pspec values could come about 
through merging, which may therefore be the most important mech- 



anism in the formation of NCC systems. It seems that there is a 
complex interplay between the gas and group properties in these 
systems. 



6 CONCLUSIONS 

We have constructed a sample of 28 galaxy groups fr om the Two- 
Dime nsional XMM-Newton Group Survey (Finoguenov et al.l2006l 
|2007|) and the RASSCALS groups of Mahdavi et al. (2005^). We 
have statistically analysed the results of a novel two-dime nsiona l 
spectral analysis techn ique applied by [ Mahdavi e t akl ( l2005h . 
iFinoguenov et alj ( l2006h and lFinoguenov et al.i Ii2007i) to high qual- 
ity XMM-Newton data. This is the largest group sample that has 
been analysed in a consistent manner with XMM-Newton data, and 
the size of the sample allowed a division on the basis of the pres- 
ence of a cool core. We find 18 groups to exhibit cool cores (CCs) 
and 10 to be classified as non cool cores (NCCs). This latter sample 
is the first of its kind in the group regime, and the analysis of these 
groups provides useful insights into the variation in the observed 
properties of these systems. We summarise our main results below: 

(i) We have measured the ratio of the central temperature drop 
(defined between the peak temperature and that at O.OI rsoo) to the 
peak temperature in the CC systems. This fractional central decline 
in T is smaller than in clusters. Within our whole group sample, 
we find no significant trend with mean temperature, however, the 
hottest system in our sample (NGC 4073, f = 1.87 keV) has an 
unusually small temperature decline, and if this group is excluded 
we find a positive correlation between CC strength and T which is 
significant at the > la level. 

(ii) We find radia lly smaller CCs in the coole st systems, con- 
firming the results o flRasmussen & PonmanI ( |2007|) . and find no CC 
systems with temperatures less than ~0.7 keV, indicating a possi- 
ble lack of CC systems at lower temperatures. 

(iii) We have investigated the relationship between temperature 
a nd entropy evaluated at 0. 1 rsoo, incorporating the cluster sample 
of'Sanderson et al.'|l2009^. Applying a BCES orthogonal regression 
(Akritas & Bershady 1996) yield s a slope of 0.7 9=fc0.06, in agree- 
ment with the Chandra sample of lSun et"al] ( 12008 ). There is a large 
amount of scatter in the group data, and the fit is mostly constrained 
by the addition of the cluster points, leading us to adopt the mod- 
ified entropy scaling (5* oc T^/^ of IPonman et al.ll2003l) for pur- 
poses of scaling group properties in this paper. 

(iv) The entropy profiles of cool groups (defined as those with 
mean temperature below the median value for our sample of 
1.035 keV) are significantly steeper than those of hotter groups. 
The slope of the entropy profiles of the Cool subsample is found 
to be 0.77±0.03, compared to 0.65±0.02 for the Warm group sub- 
sample. 

(v) Comparing the entropy profiles of CC and NCC groups, we 
find there to be more scatter in the entropy profiles for the latter, 
which is consistent with the larger scatter shown in their tempera- 
ture profiles. The central entropy, within ~0. 1 7-500, is also flatter 
for the NCC systems, whereas in CCs, the entropy profiles decrease 
steadily into the centre of the system. There is no evidence for a 
central entropy pedestal, as required for clusters dOonahue et al.l 
|2006|) . in the CC groups. 

(vi) Examining the pressure and density profiles, we find a 
marked difference in the pressure profiles of Cool and Warm 
groups, with the scaled pressure profile of the former rising to 
higher central values than in warmer systems. We attribute this dif- 
ference to the increased dominance of the stellar mass of the bright- 
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est group galaxy in the centres of the cooler systems. Comparing 
CC and NCC systems, we see slightly steeper density and pressure 
profiles within the cores of CC systems, which is attributable to 
their central temperature decline. Outside the typical radius of the 
CC (~0.1r5oo), we find the two subsamples to have similar gas 
profiles. 

(vii) We have investigated the entropy distribution of the 
groups by plotting the observed gas mass in bins of scaled 
entropy. The total gas mass within 0.5 rsoo is compara- 
ble for the CC (4.7±0.5x lO^M© keV"^) and NCC groups 
(4.1±0.5x lO^^M© keV~^). Comparing the observed entropy dis- 
tributions to the theoretical expectation if non-gravitational pro- 
cesses are ignored, and find that the observed distributions gen- 
erally reach much higher entropies within O.Srsoo. In trying to 
reconcile the differences between the observed and theoretical his- 
tograms, we find that simple modifications (shift, truncation or ra- 
diative cooling) are not sufficient to bring the theoretical distri- 
bution in line with that observed. A 'shift -I- cool' model which 
aims to match the high entropy behaviour of the two distributions 
performs well when the entropy shift is small, but in many cases 
the large entropy boost suppresses cooling, resulting in an entropy 
profile which lacks the low entropy gas required by the observa- 
tions. This suggests a process whereby the entropy shift is larger 
for higher entropy gas. Potential mechanisms for achieving this are 
either through entropy amplification or episodic heating. We note 
that the former is likely to be more energy efficient. 

(viii) We define the 'feedback impact' of a group, and find that 
systems with the highest feedback impact reach higher scaled en- 
tropies within 0.5 rsoo, regardless of whether or not they have CCs. 
We find higher metallicity in the central regions of the lowest feed- 
back systems, but over the majority of the radial range, the metal- 
licity profiles of the low and high feedback systems is comparable. 
If low entropy metal-rich gas in the highest feedback systems has 
been boosted in entropy, it has been pushed outside 0.5 r^oo- There 
is no evidence for an increase in the metal content with the level 
of feedback impact, which leads us to favour AGN, probably act- 
ing before group assembly, as the dominant source of feedback, 
rather than supernova explosions. We test for bias in our defini- 
tion of feedback impact by scaling by the total mass of the system 
(approximated as oc T^/^) rather than the gas mass. Changing our 
definition of feedback impact in this way does not affect our results. 

(ix) Fitting the a - T relation to the group sample yields a steep 
slope of 0.9±0.3. Six of the NCC groups lie above the l3spec=l 
line, which seems inconsistent with models which invoke feedback 
to eliminate CCs in these systems. We suggest that group-group 
mergers are more likely responsible for the elimination of CCs. 
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